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Introduction 


X-ray  mammography  currently  serves  as  the  primary  screening  method  for  breast  cancer.  This  research 
seeks  to  improve  the  efficacy  of  mammography  by  optimizing  the  entire  imaging  chain  for  the  detection  of 
breast  masses  and  calcifications.  Our  research  can  be  decomposed  into  two  distinct  areas.  The  first 
area  involves  a  thorough  understanding  of  the  physical  characteristics  of  the  imaging  chain.  This 
includes  experimental  measurements  of  the  resolution  and  noise  of  x-ray  detectors  and  display  devices. 

In  addition,  this  requires  modeling  of  physical  processes,  such  as  scatter,  that  impact  image  quality.  The 
second  area  applies  this  physical  data  to  simulation  studies,  allowing  one  to  explore  a  variety  of 
parameters  in  the  imaging  chain.  Finally,  the  simulation  procedures  create  images  that  simulate  the 
resolution  and  noise  of  various  imaging  configurations.  A  combination  of  observer  models  and  human 
observers  view  these  simulated  images  in  order  to  discover  the  impact  of  these  resolution  and  noise 
settings  on  the  detection  of  breast  masses  and  calcifications. 

Body 

This  section  reviews  the  progress  of  the  research  in  addressing  the  approved  statement  of  work.  We 
have  included  only  those  parts  of  the  statement  of  work  addressed  by  our  work  over  the  past  year. 
Sections  not  included  in  this  report  are  part  of  our  future  work. 

Task  1:  Create  a  simulation  procedure  for  the  anatomical  background  of  mammographic  images 

1. 1  Acquire  normal  mammograms  obtained  on  digital  systems  for  analysis 

Working  with  colleagues  from  Emory  University,  we  obtained  984  images  acquired  on  an  indirect 
flat-panel  detector. 

1.2  Categorize  the  images  into  the  four  types  of  breast  composition,  as  identified  by  the  BIRADS 
system. 

1.3  Analyze  the  geometrical  features  of  these  breasts  and  characterize  them  with  a  fixed  number  of 
scalar  parameters,  such  as  size. 

These  two  steps  were  included  as  they  would  aid  in  the  creation  on  a  routine  to  simulate 
mammographic  backgrounds.  As  part  of  the  research  for  anatomical  simulation,  we  searched  the 
literature  for  previous  research  on  mammographic  background  simulation.  We  discovered  and 
implemented  the  methods  of  Bochud,  ef  a/to  emulate  mammographic  backgrounds  by  creating 
clustered  lumpy  backgrounds.1  These  simulated  backgrounds  appeared  similar  to  real 
mammographic  backgrounds,  but  did  not  capture  all  of  the  complexity  of  real  anatomy. 

Therefore,  we  decided  to  use  the  mammographic  data  set  obtained  in  1.1  for  our  subsequent 
simulation  experiments. 

1.4  Obtain  mammograms  from  the  Digital  Database  for  Screening  Mammography  (DDSM)  to  analyze 
lesion  characteristics 

We  selected  images  from  the  DDSM  that  contained  oval  circumscribed,  oval  obscured,  irregular 
ill-defined,  and  irregular  spiculated  masses.  In  addition,  we  selected  images  that  contained  fine 
linear  branching  and  pleomorphic  calcifications.  We  segmented  these  mammograms  into  regions 
of  2.56  cm  x  2.56  cm  centered  on  the  mass  or  calcification. 

1.5  Analyze  the  features  of  specific  lesion  types 

We  analyzed  the  images  obtained  from  the  DDSM  to  create  a  model  of  the  radiographic 
appearance  of  breast  lesions.  This  model  was  described  in  our  2004  IWDM  proceedings  article 
in  Appendix  III. 

1.6  Create  a  program  that  can  create  images  with  breast  anatomy  and  breast  lesions  that  allows  for 
user  Input  of  specific  scalar  parameters,  such  as  size. 

We  created  a  program  that  allowed  one  to  insert  simulated  masses  and  calcifications  into  normal 
anatomical  backgrounds.  The  details  of  this  program  were  disclosed  in  our  2004  IWDM 
proceedings  article  in  Appendix  III.  In  addition,  this  program  served  a  crucial  role  in  our  recent 
research  on  the  impact  of  display  image  quality,  as  disclosed  in  Appendix  II. 

As  noted  in  1 .2  and  1 .3,  the  simulated  mammographic  backgrounds  lacked  the  complexity  of  real 
backgrounds.  We  therefore  used  actual  digital  mammograms  for  our  simulation  experiments. 
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1. 7  Establish  mapping  technique  to  determine  grayscale  values  of  image  using  sigmoid  curve 
transformation. 

To  conduct  this  task,  an  experienced  mammographer  reviewed  the  digital  images  obtained  in  1 .1 
and  window  and  leveled  each  mammogram  to  produce  a  clinically  relevant  appearance.  We 
recorded  the  parameters  for  each  image  and  fit  a  sigmoid  curve  to  each  window  and  level 
function.  We  then  applied  the  appropriate  transformation  to  each  image  in  order  to  simulate  the 
correct  clinical  appearance.  This  stage  was  disclosed  in  our  2005  SPIE  proceedings  article  in 
Appendix  II. 

Task  3:  Create  an  empirical  model  that  relates  the  resolution  and  noise  of  a  digital  mammographic 
system  to  the  detectability  of  breast  lesions. 

3. 1  Compile  a  list  of  MTFs  and  NPS  for  commercial  radiographic  systems,  including  image 
processing  algorithms  and  displays. 

We  have  begun  compiling  a  list  of  MTFs  and  NPSs  for  several  mammographic  systems.  In 
addition,  we  have  conducted  some  studies  to  directly  measure  the  physical  characteristics  of 
mammographic  systems.  Please  refer  to  Appendix  I  for  one  study  where  we  measured  the 
performance  of  a  clinical  prototype  digital  mammographic  system. 

Task  5:  Utilize  the  empirical  model  to  examine  the  effect  of  dose  on  the  detection  of  microcalcifications 
and  masses  and  determine  the  minimum  allowable  dose  level  for  "safe”  mammographic  imaging. 

5. 1  Determine  the  relationship  between  dose  and  noise  amplitude  for  the  three  specific  digital 
mammographic  systems  through  published  measurements. 

We  determined  the  magnitude  of  the  signal  to  noise  ratio  for  a  given  dose  by  the  equation: 
SNRAclJ=DQE(0)-SNRlM2  (1) 

where  SNR,deat  was  computed  using  a  program  by  Boone  to  generate  x-ray  spectra2  and  DQE( 0) 
was  determined  from  published  measurements.  This  signal  to  noise  ratio  was  mapped  to  a 
graylevel  variance  using  the  exposure-pixel  value  relationship  for  the  detector. 

5.2  Determine  the  effect  of  scatter  utilizing  previously  published  models. 

We  determined  the  magnitude  of  scatter  by  using  previously  published  data  by  Boone.3  Our 
group  measured  the  magnitude  of  scatter  reduction  accomplished  by  the  antiscatter  grid.  The 
scatter  to  primary  ratios  were  then  discounted  by  the  scatter  reduction  from  the  grid.  The  effect  of 
scatter  was  incorporated  by  reducing  the  contrast  of  our  simulated  lesions  by  the  magnitude  of 
the  scattered  radiation. 

5.3  Using  the  previously  developed  empirical  model  to  analyze  the  effect  of  dose  on  the  detectability 
of  masses  and  microcalcifications. 

We  have  generated  four  image  sets  using  the  mammographic  data  obtained  in  1 .1 .  The  first  set 
was  obtained  at  full  dose  and  the  next  three  sets  have  added  noise  to  simulate  half,  quarter,  and 
eighth  dose,  respectively.  Two  different  observer  models  analyzed  these  image  sets,  a  visual 
discrimination  model  and  a  non-prewhitening  matched  filter  with  eye  filter  model.  We  are 
currently  examining  the  results  from  these  observer  models. 

Task  6:  Apply  the  empirical  model  to  ascertain  the  effect  of  a  specific  image  processing  algorithm, 
unsharp  masking,  on  lesion  detection  and  optimize  its  utilization. 

6. 1  Examine  the  clinical  parameters  used  for  unsharp  masking. 

Several  types  of  unsharp  masking  are  used  in  clinical  practice.  We  implemented  the  most  basic 
type  of  multiscale  processing,  consisting  of  an  unsharp  masking  stage  and  a  contrast  equalization 
stage.  The  form  of  this  processing  and  parameters  were  determined  from  previously  published 
methods.4 

Task  7:  Employ  the  model  to  examine  the  influence  of  two  specific  display  characteristics,  display 
magnification  and  display  resolution,  on  lesion  detection  and  thus  develop  guidelines  for  optimized 
viewing  of  digital  mammograms. 

7.2  Determine  the  resolution  and  noise  for  four  display  devices,  three  common  Cathode  Ray  Tube 
(CRT)  devices  and  one  Liquid  Crystal  Display  (LCD)  device. 
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Using  a  high-quality  CCD  camera,  we  measured  the  resolution  and  noise  of  two  CRT  displays 
and  three  LCD  devices.  As  LCDs  are  becoming  increasingly  common  in  clinical  systems,  we 
decided  to  include  more  of  a  focus  on  LCD  displays  than  we  proposed  in  the  statement  of  work. 
We  have  written  a  manuscript  describing  our  measurement  protocol  and  experimental  results  and 
submitted  it  to  Medical  Physics  in  mid-March  2005. 

7.3  Fit  the  resolution  and  noise  properties  of  the  combined  display  and  detector  system  using  the 
generalized  curve-fitting  algorithm. 

After  obtaining  the  resolution  and  noise  characteristics,  we  fit  each  of  them  with  a  multi-parameter 
exponential  function.  This  provided  us  with  a  functional  form  for  the  resolution  and  noise  data. 

7.4  Input  the  above  into  the  empirical  model  in  order  to  develop  guidelines  for  optimized  display  of 
mammographic  images. 

Instead  of  an  empirical  model,  we  used  observer  models  to  examine  the  impact  of  different 
display  resolutions  on  the  detection  of  masses  and  calcifications.  Please  refer  to  our  2005  SPIE 
proceedings  article  in  Appendix  II  for  the  full  details  of  this  project. 
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Key  Research  Accomplishments 

•  Acquired  a  large  data  set  of  normal  digital  mammograms. 

•  Developed  model  for  radiographic  appearance  of  breast  masses  and  calcifications  and 
implemented  lesion  simulation  program. 

•  Measured  resolution  and  noise  of  five  medical  displays,  representing  both  CRT  and  LCD  devices. 

•  Created  large  image  set  with  resolution  properties  emulating  those  of  commercial  medical 
displays. 

•  Created  large  image  set  that  simulated  the  effects  of  reduced  dose. 

•  Implemented  observer  model  for  examining  image  sets,  based  on  a  non-prewhitening  matched 
filter  model  with  eye  filter. 

•  Measured  physical  characteristics  of  clinical  prototype  mammographic  system. 

•  Researched  image  processing  techniques  and  wrote  program  based  on  a  common  image 
processing  algorithm. 
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Reportable  Outcomes 

CONFERENCE  PRESENTATIONS 

R.S.  Saunders  Jr.,  E.  Samei,  J.  Johnson,  and  J.  Baker,  “Effect  of  Display  Resolution 

on  the  Detection  of  Mammographic  Lesions,”  SPIE  Medical  Imaging  2005:  Image  Perception,  Observer 

Performance,  and  Technology  Assessment  (2005). 

R.  S.  Saunders  Jr.,  E.  Samei,  J.  Baker,  “Simulation  of  Mammographic  Lesions,”  Radiological  Society  of 
North  America  Annual  Meeting  (2004). 

R.  S.  Saunders  Jr,  A.  Farshchi,  E.  Samei,  “Measurement  of  Display  Resolution  for  Commercial  Medical 
Displays,”  Radiological  Society  of  North  America  Annual  Meeting  (2004). 

R.  S.  Saunders  Jr.,  E.  Samei,  J.  Y.  Lo  and  J.  L.  Jesneck,  "Physical  Characterization  of  a  Selenium-based 
Full  Field  Digital  Mammography  Detector,"  7th  International  Workshop  on  Digital  Mammography  (2004). 

R.  S.  Saunders  Jr.,  E.  Samei,  and  J.  Baker,  “Simulation  of  Breast  Lesions,”  7th  International  Workshop  on 
Digital  Mammography  (2004). 

REFEREED  JOURNAL  ARTICLES 

R.S.  Saunders  Jr.,  E.  Samei,  J.L.  Jesneck,  and  J.Y.  Lo,  “Physical  characterization  of  a  prototype 
selenium-based  full  field  digital  mammography  detector,”  Med.  Phys.  32,  588-599  (2005). 
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Conclusions 


This  year,  we  have  developed  several  tools  that  allow  us  to  explore  the  impact  of  resolution  and  noise  on 
the  detection  of  mammographic  lesions.  First,  we  have  developed  a  routine  that  inserts  simulated 
masses  or  calcifications  into  a  normal  mammographic  background.  Second,  we  have  measured  the 
physical  properties  of  a  clinical  mammographic  detector  and  several  medical  display  devices.  Third,  we 
have  acquired  a  large  data  set  of  mammographic  images.  Applying  these  tools,  we  have  explored  two 
questions  this  year.  The  first  question  explored  the  impact  of  display  resolution  on  the  detection  of  breast 
masses  and  calcifications.  The  second  question  explored  the  effect  of  reduced  dose  on  the  detection  of 
breast  lesions.  Both  of  these  questions  have  immediate  impact  on  clinical  care,  as  they  will  determine 
which  medical  displays  are  appropriate  for  reading  mammograms  and  whether  women  may  be  imaged 
using  a  lower  dose.  Future  work  will  include  examining  the  effect  of  image  processing  and  the  relative 
merit  of  different  mammographic  detectors  on  the  detection  of  breast  lesions. 
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(Received  23  July  2004;  revised  15  November  2004;  accepted  for  publication  14  December  2004; 
published  3  February  2005) 

The  purpose  of  this  study  was  to  measure  experimentally  the  physical  performance  of  a  prototype 
mammographic  imager  based  on  a  direct  detection,  flat-panel  array  design  employing  an  amorphous 
selenium  converter  with  70  pm  pixels.  The  system  was  characterized  for  two  different  anode  types, 
a  molybdenum  target  with  molybdenum  filtration  (Mo/ Mo)  and  a  tungsten  target  with  rhodium 
filtration  (W/Rh),  at  two  different  energies,  28  and  35  kVp,  with  approximately  2  mm  added 
aluminum  filtration.  To  measure  the  resolution,  the  presampled  modulation  transfer  function  (MTF) 
was  measured  using  an  edge  method.  The  normalized  noise  power  spectrum  (NNPS)  was  measured 
by  two-dimensional  Fourier  analysis  of  uniformly  exposed  mammograms.  The  detective  quantum 
efficiencies  (DQEs)  were  computed  from  the  MTFs,  the  NNPSs,  and  theoretical  ideal  signal  to 
noise  ratios.  The  MTF  was  found  to  be  close  to  its  ideal  limit  and  reached  0.2  at  1 1.8  mm-1  and  0.1 
at  14.1  mm'1  for  images  acquired  at  an  RQA-M2  technique  (Mo/Mo  anode,  28  kVp,  2  mm  Al). 

Using  a  tungsten  technique  (MW2;  W/Rh  anode,  28  kVp,  2  mm  Al),  the  MTF  went  to  0.2  at 
11.2  mm'1  and  to  0.1  at  13.3  mm'1.  The  DQE  reached  a  maximum  value  of  54%  at  1.35  mm'1  for 
the  RQA-M2  technique  at  1.6  yu-C/kg  and  achieved  a  peak  value  of  64%  at  1.75  mm'1  for  the 
tungsten  technique  (MW2)  at  1.9  /xC/kg.  Nevertheless,  the  DQE  showed  strong  exposure  and 
frequency  dependencies.  The  results  indicated  that  the  detector  offered  high  MTFs  and  DQEs,  but 
structured  noise  effects  may  require  improved  calibration  before  clinical  implementation.  ©  2005 
American  Association  of  Physicists  in  Medicine.  [DOI:  10.1118/1.1855033] 
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I.  INTRODUCTION 

Breast  cancer  remains  the  second  leading  cause  of  cancer 
death  for  women  in  the  United  States.  The  American  Cancer 
Society  (ACS)  estimates  that  in  2004,  215  990  new  cases  of 
invasive  breast  cancer  will  be  diagnosed  and  40  110  women 
will  die  from  the  disease  in  the  United  States.1  Early  detec¬ 
tion  of  this  disease  holds  the  key  for  survival,  as  more  treat¬ 
ment  options  exist  for  early  stage  cancers  and  treatments 
tend  to  be  more  successful  at  this  stage.  X-ray  mammogra¬ 
phy  continues  to  be  widely  regarded  as  the  most  effective 
early-detection  screening  tool  available  today.2,3  X-ray  mam¬ 
mography  places  severe  demands,  however,  on  an  imaging 
system.  A  system  must  capture  small,  low  contrast  anatomi¬ 
cal  details,  as  the  early  signs  of  cancer  are  often  very  subtle. 
While  mammography  has  experienced  notable  advancements 


in  recent  years,  further  improvement  is  required  as  up  to  22% 
of  cancers  are  missed  at  the  initial  screening.4 

Full  Field  Digital  Mammography  (FFDM)  offers  the 
promise  of  improving  mammographic  image  quality  and 
therefore  increasing  the  utility  of  this  screening  procedure.5'7 
As  images  are  stored  in  a  digital  format,  a  radiologist  can 
view  the  images  at  any  workstation  or  many  clinicians  can 
have  simultaneous  access  to  the  images.  The  use  of  image 
processing  algorithms  enhances  various  features  in  the  im¬ 
age.  In  addition,  these  systems  have  the  potential  to  improve 
mammographic  imaging  by  separating  each  stage  of  the  im¬ 
aging  chain,  from  detection  to  image  processing  to  display, 
allowing  each  step  to  be  independently  optimized. 

The  current  state  of  the  art  in  digital  mammography  is 
solid-state  flat-panel  detectors.8  Flat-panel  detectors  can  be 
subdivided  into  two  categories,  direct  and  indirect,  named 
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for  the  mechanisms  used  to  detect  x-rays.9,10  In  direct  detec¬ 
tors,  a  photoconductive  layer  absorbs  an  incoming  x-ray  pho¬ 
ton  and  converts  it  to  electric  charge.  A  voltage  applied 
across  the  photoconductor  then  draws  the  charges  toward  the 
pixel  electrodes.11'12  In  contrast,  indirect  detectors  utilize  a 
scintillation  layer  that  converts  the  x-ray  photon  into  visible 
light  photons,  which  are  subsequently  absorbed  by  photosen¬ 
sitive  elements.13,14  Because  of  the  different  physical  mecha¬ 
nisms  used  to  detect  photons,  the  image  quality  characteris¬ 
tics  of  these  detectors  differ  substantially.  Several  prior 
studies  have  substantiated  some  of  these  differences.9,1  ,15,16 
In  addition,  two  previous  studies  have  examined  limited  as¬ 
pects  of  image  quality  for  selected  mammographic  detectors 
using  amorphous  selenium.17,18 

The  main  purpose  of  this  work  was  to  comprehensively 
evaluate  the  physical  image  quality  characteristics  of  an  early 
prototype  mammographic  detector  based  on  a  direct  detec¬ 
tion  flat-panel  array  design  that  employed  an  amorphous  se¬ 
lenium  converter.  Three  key  metrics  of  image  quality  were 
evaluated  for  several  radiographic  techniques,  the  modula¬ 
tion  transfer  function  (MTF),  normalized  noise  power  spec¬ 
trum  (NNPS),  and  detective  quantum  efficiency  (DQE), 
which  described  the  resolution,  noise,  and  signal  to  noise 
performance  of  the  detector,  respectively. 19-24  As  previous 
research  has  shown  that  selenium  detectors  can  exhibit  im¬ 
age  lag  and  ghosting,25  this  research  also  examined  the  lag 
performance  of  the  detector. 

A  secondary  objective  of  this  research  was  to  consider 
new  beam  qualities  for  digital  mammography.  Traditionally, 
screen-film  mammography  was  performed  using  a  beam 
from  a  molybdenum  target  with  molybdenum  filtration.26 
This  beam  quality  might  not  be  optimal  for  digital  mammog¬ 
raphy,  however,  given  the  different  energy  sensitivities  and 
greater  dynamic  range  of  digital  detectors.  Several  research¬ 
ers  had  suggested  that  other  beam  qualities  could  allow  for 
better  image  quality  for  digital  mammography.  ’  There¬ 
fore,  the  study  examined  the  image  quality  characteristics  for 
two  different  anode  types,  a  molybdenum  target  with  molyb¬ 
denum  filtration  and  a  tungsten  target  with  rhodium  filtration, 
and  for  two  different  energies,  28  kVp  and  35  kVp,  with 
added  aluminum  filtration. 

II.  METHODS  AND  MATERIALS 
A.  Detector  description 

The  detector  investigated  in  this  study  was  an  early  pro¬ 
totype  mammographic  imager  based  on  a  direct  detection 
flat-panel  array  design  that  employed  an  amorphous  sele¬ 
nium  converter  (Mammomat  NovationDR;  Siemens  Medical 
Solutions;  Erlangen,  Germany).  The  detector  utilized  a 
250  pm  amorphous  selenium  photoconductive  layer  coupled 
to  a  matrix  of  pixels,  each  with  a  storage  capacitor  and  amor- 
phous  silicon  switching  transistor.  The  active  detector  area 
was  23.3  cm  X  28.7  cm  consisting  of  3328X4096  square 
pixels.  Each  pixel  was  placed  with  a  70  pm  pixel  pitch  and 
offered  a  fill  factor  of  greater  than  90%.  This  product  has 
since  received  FDA  approval. 


Fig.  1.  Coordinate  system  for  physical  measurements.  These  axes  are  la¬ 
beled  by  the  anatomy  imaged  in  the  craniocaudal  view. 

Prior  to  evaluation,  the  standard  antiscatter  grid  and  com¬ 
pression  paddle  were  removed  from  the  system.  For  most 
measurements,  the  standard  detector  cover  was  placed  on  the 
system.  For  the  MTF  measurements,  the  detector  cover  was 
removed  so  that  an  edge  device  could  be  placed  as  close  as 
possible  to  the  active  selenium  layer  to  minimize  focal  spot 
blur. 

The  coordinate  system  used  to  describe  the  system,  as 
shown  in  Fig.  1,  referred  to  the  anatomical  features  as 
viewed  on  a  craniocaudal  view.  There  were  two  main  axes, 
the  chest  wall-nipple  (CN)  axis  as  well  as  the  left-right  (LR) 
axis.  By  examining  the  system  performance  along  these  two 
orthogonal  axes,  one  was  able  to  identify  any  asymmetries. 

B.  Image  acquisition 

A  high-frequency,  multiphase  x-ray  generator  (Mammo¬ 
mat  NovationDR),  for  which  the  high  voltage  accuracy  was 
verified  to  be  within  ±5%,  served  as  the  x-ray  source  for  the 
system.  The  anode  was  operated  with  a  large  focal  spot  of 
0.3  mm  (IEC),  nominal,  for  all  image  acquisitions.  No  post¬ 
processing  was  applied  to  the  images.  All  images  were  trans¬ 
ferred  to  a  research  computer  as  14-bit,  raw  data  for  analysis. 

Prior  to  image  acquisition,  the  detector  underwent  routine 
detector  calibration  to  correct  for  dead  pixels  and  gain  non¬ 
uniformities.  The  process  formed  a  dead  pixel  map  by  de¬ 
tecting  inactive  pixels  in  a  flat-field  image  acquired  at  28 
kVp  with  a  4  cm  PMMA  slab  in  the  beam.  A  gain  map  was 
similarly  computed  from  the  average  of  eight  flat-field  im¬ 
ages  also  acquired  at  28  kVp  with  a  4  cm  PMMA  slab  in  the 
beam.  As  no  images  in  this  research  utilized  an  antiscatter 
grid,  the  calibration  was  performed  without  a  grid  in  place. 
The  system  corrected  all  subsequently  acquired  images  using 
the  gain  and  dead  pixel  maps. 

For  all  image  acquisitions,  the  exposure  to  the  detector 
was  measured  using  a  calibrated  ionization  chamber  (1515 
x-ray  monitor  with  10X5-6M  dedicated  mammography  ion¬ 
ization  chamber,  Radcal  Corporation,  Monrovia,  CA)  placed 
at  48  cm  from  the  focal  spot.  As  reported  in  previous  studies, 
this  ionization  chamber  had  little  energy  dependence  over 
mammographic  energies.31  Manufacturer  specifications  note 
that  the  calibration  accuracy  of  the  chamber  was  ±4%  (at  20 
kVp,  0.26  mm  A1  HVL)  with  ±5%  energy  dependence  in  the 
10  keV  to  40  keV  range.  The  exposures  incident  on  the  de¬ 
tector,  located  at  65  cm  distance  from  the  focal  spot,  were 
estimated  from  the  measured  exposure  values  using  the 
inverse-square  law. 
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Table  I.  Beam  qualities  used  for  physical  characterization  of  the  detector. 
The  aluminum  used  for  the  added  filtration  had  s=99%  purity. 


Name 

Anode  target 

Anode 

filtration 

kVp 

Added 
filtration 
(mm  Al) 

Half-value 
layer 
(mm  Al) 

RQA-M2 

Molybdenum 

Molybdenum 
(30  urn) 

28 

2 

0.6 

RQA-M4 

Molybdenum 

Molybdenum 
(30  |rm) 

35 

1.8 

0.68 

MW2 

Tungsten 

Rhodium 
(50  pm) 

28 

2 

0.79 

MW4 

Tungsten 

Rhodium 
(50  pm) 

35 

2 

0.92 

mine  the  presampled  modulation  transfer  function  (MTF).  In 
summary,  first  a  double  Radon  transformation  determined  the 
angle  of  the  edge  transition  with  0.01°  accuracy.  The  edge 
spread  function  (ESF)  was  computed  by  projecting  the  image 
data  along  lines  parallel  to  the  edge  transition  using  bin  sizes 
of  0.1  pixels.  To  minimize  noise,  the  ESF  was  smoothed 
using  a  modest  fourth-order  moving  polynomial  fit  and  dif¬ 
ferentiated  to  form  the  line  spread  function  (LSF).  A  Han¬ 
ning  window  with  10  mm  width  was  then  applied  to  the  LSF 
to  force  the  tails  of  the  LSF  to  zero.  Finally,  the  presampled 
MTF  was  computed  as  the  normalized  Fast  Fourier  Trans¬ 
form  of  the  LSF. 


Four  different  beam  qualities  were  utilized  for  the  image 
quality  measurements,  as  outlined  in  Table  I.  Two  molybde¬ 
num  techniques,  RQA-M2  and  RQA-M4,  were  chosen  from 
the  International  Electrotechnical  Commission  (IEC)  stan¬ 
dard  6 1267-2. 32  The  standard  specified  the  anode  type,  anode 
filtration,  kVp,  and  half-value  layer  for  each  beam  quality. 
Aluminum  filtration  was  then  placed  in  the  beam  to  produce 
the  desired  half-value  layer.  The  IEC  standard  did  not  in¬ 
clude  corresponding  tungsten  techniques  for  mammographic 
applications.  To  facilitate  meaningful  comparisons  between 
detector  systems,  two  additional  tungsten  techniques,  MW2 
and  MW4,  were  used  that  had  similar  characteristics  to  those 
of  the  molybdenum  techniques,  RQA-M2  and  RQA-M4.  The 
half-value  layers  for  these  four  beam  qualities  were  mea¬ 
sured  using  a  narrow  geometry  and  added  aluminum  filtra¬ 
tion  in  0.1  mm  increments  around  the  estimated  half-value 
layer  thickness.  The  half-value  layer  thicknesses  were  then 
estimated  from  logarithmic  interpolation  of  the  measured  ex¬ 
posure  values.32 

C.  Linearity 

Linearity  was  determined  by  exposing  the  detector  to  a 
wide  range  of  uniform  x-ray  exposures  for  each  of  the  four 
radiographic  techniques  described  above.  The  average  pixel 
values  were  computed  from  a  14.3  cmX  14.3  cm  region  lo¬ 
cated  near  the  chest  wall  section  of  the  detector.  From  this, 
the  relationships  between  mean  pixel  value  and  exposure 
were  ascertained  for  each  technique. 

D.  Modulation  transfer  function 

An  edge  method,  reported  in  previous  publi¬ 
cations, 9,I0,33“36  was  used  to  measure  the  presampled  MTF.  A 
0. 1  mm  Pt-Ir  edge  was  placed  in  contact  with  the  detector  at 
1  cm  distance  from  the  chest  wall  edge  of  the  detector.  The 
device  was  oriented  at  a  3°-6°  angle  with  respect  to  the  pixel 
array.  Edge  images  were  then  acquired  at  each  of  the  four 
radiographic  techniques  at  relatively  high  exposure  values  of 
16.2  /xC/kg  (62.6  mR),  15.3  fiC/kg  (59.2  mR),  9.52  /x,C/kg 
(36.9  mR),  and  9.75  fiC/kg  (37.8  mR)  for  RQA-M2, 
RQA-M4,  MW2,  and  MW4  techniques,  respectively. 

A  previously  reported  routine  analyzed  the  edge  images 
in  a  region  around  the  edge  (21.2  mm  X  35.8  mm)  to  deter- 


E.  Normalized  noise  power  spectrum 


To  characterize  the  system  noise,  images  were  acquired  of 
uniform  beams  of  radiation  for  the  different  techniques, 
while  the  exposure  was  simultaneously  measured  with  an 
ionization  chamber.  The  NNPS  was  then  computed  from 
these  flat-field  images  using  previously  published 
methods.36-38  A  large  region  near  the  chest  wall  side  of  the 
detector,  excluding  the  edges  of  the  image,  was  used  for 
analysis.  This  region  was  segmented  into  256  sequential  re¬ 
gions  of  interest  (ROIs)  of  128  X  128  pixels.  A  two- 
dimensional  polynomial  surface  was  subtracted  from  each 
region  of  interest  (ROl)  to  minimize  background  trending 
and  a  Hamming  window  was  applied  to  each  ROI  so  that  the 
edges  of  the  ROI  went  to  zero.  To  account  for  intensity  varia¬ 
tions  in  the  image,  each  ROI  was  then  scaled  by  the  ratiq  of 
its  mean  to  the  mean  pixel  value  of  the  ROI  in  the  top-left- 
hand  corner  of  the  image.  Each  ROI  was  transformed  by  a 
two-dimensional  FFT  and  the  absolute  magnitude  squared  of 
each  FFT  was  averaged  together  to  obtain  the  NNPS.  This 
procedure  could  be  summarized  in  the  following  equation:37 


NNPS(w,i>) 


M  ( 

dA  y  I  (ROI,) 
M-N2~[  [ (ROI|> 


FFT 


1 

.(ROI,) 


(ROI; 


(1) 


where  dA  represented  the  pixel  area,  M  described  the  number 
of  regions  of  interest  in  which  the  image  was  segmented,  N 
corresponded  to  the  number  of  pixels  along  one  edge  of  an 
ROI,  ROI,-  referred  to  a  particular  region  of  interest  within 
the  flat  field  image,  ROI;  corresponded  to  the  ROI  in  the 
top-left  comer  of  the  image,  and  (ROI,)  was  the  mean  of 
ROI,-.  To  summarize  this  two-dimensional  information  in 
one-dimensional  form,  horizontal  and  vertical  traces  were 
obtained  by  averaging  together  the  central  frequency  bands 
(the  central  axis  and  ±5  frequency  lines).  Radial  traces  were 
also  obtained  by  radial  averaging. 

The  magnitude  of  the  NNPS  could  be  related  to  the  image 
variance  using  Parseval’s  Theorem39  and  applying  ergodic 
assumptions.  This  allowed  the  replacement  of  (ROI,)  by  (/), 
the  mean  of  the  entire  image,  and  the  mean  variance  of  the 
ROIs  became  the  variance  of  the  image,  a2.  One  could  then 
show  that 
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Table  II.  Ideal  SNR2/mR  values  calculated  for  an  energy-integrating  detec¬ 
tor.  The  beams  were  modeled  with  the  specified  intrinsic  filtrations  as  well 
as  the  experimentally  measured  half-value  layer. 


Name 

Anode  target 

Anode  filtration 

9,deai(mm  2  mR-1) 

RQA-M2 

Molybdenum 

Molybdenum  (30  pm) 

46052 

RQA-M4 

Molybdenum 

Molybdenum  (30  pm) 

52542 

MW2 

Tungsten 

Rhodium  (50  pm) 

54773 

MW4 

Tungsten 

Rhodium  (50  pm) 

67781 

^  ,  N2dA  , 

2NNPS(«,i>)  =  — -o2.  (2) 

u,v  W 

For  a  linear,  quantum-limited  detector,  (I)  and  o2  are  propor¬ 
tional  to  the  exposure,  E,  which  would  make  the  product  of 
the  NNPS  and  exposure  independent  of  exposure.  The  prod¬ 
uct  of  NNPS  and  exposure  was  then  used  as  a  way  to  assess 
how  well  the  detector  approximated  a  quantum-limited  de¬ 
tector. 

A  second  examination  of  system  noise  utilized  a  back¬ 
ground  subtraction  method,  which  isolated  the  quantum 
noise  components  of  total  system  noise.39  42  An  average  im¬ 
age  was  created  from  ten  repeated  images  acquired  with  the 
RQA-M2  technique  at  125  mAs.  The  average  image  was 
then  subtracted  from  one  of  the  individual  images  to  form  a 
“background-free”  image.  The  NNPS  was  then  computed 
from  the  “background-free”  image.  To  correct  for  the  change 
in  image  variance  caused  by  the  averaging  technique,  the 
NNPS  for  the  “background-free”  image  was  multiplied  by 
N/(N- 1),  where  N  equaled  10,  the  number  of  images  used 
to  create  the  average  image.42 


F.  Detective  quantum  efficiency 


The  measured  MTF  and  NNPS  were  combined  to  deter¬ 
mine  the  Detective  Quantum  Efficiency  (DQE)  as 


DQE(w)  = 


MTF2(») 

^ Ideal  ’  E  ■  NNPS(m)  ’ 


(3) 


where  g]dcal  described  the  ideal  signal  to  noise  (SNR)  ratio 
squared  per  unit  exposure  for  an  energy-integrating  detector, 
and  E  represented  the  exposure  value  at  the  detector.39,43  An 
x-ray  simulation  program  (xSpect,  Henry  Ford  Health  Sys¬ 
tem)  was  used  to  calculate  the  gIdc;d  using  a  semiempirical 
model  for  the  x-ray  spectra44  and  the  attenuation  properties 
of  the  material.33  The  q  values  are  reported  in  Table  II. 


G.  Image  lag  measurement 

The  magnitude  of  the  multiplicative  image  lag  was  char¬ 
acterized  using  the  procedure  described  in  IEC  standard 
62220-1. 45  First,  an  image  was  acquired  of  a  uniform  radia¬ 
tion  field  at  a  given  exposure  at  time  A  second  image  was 
then  acquired  of  an  edge  device  at  the  same  exposure  level  at 
time  t2.  After  a  specified  delay  time  t,  a  third  image  was 
acquired  of  a  uniform  radiation  field  at  time  t3.  This  proce¬ 
dure  then  measured  the  residual  signal  from  the  edge  device 
in  the  later  image. 


Image  I  Image  2  Image  3 


£-;  .  .  VX  svvj  .  .  '• 

Wait 

|»  □ 

j  m  ^ 

Time  z 

m  □ 

■  i 

E3  ROI 1 
M  ROI  2 


Edge  Test 
Device 


Fig.  2.  Illustration  of  the  lag  measurement  procedure  as  described  in  the 
IEC  standard  62220-1. 


To  determine  the  magnitude  of  the  residual  signal,  the 
image  data  were  examined  for  two  regions  within  all  three 
images,  as  shown  graphically  in  Fig.  2.  An  ROI,  ROIb  was 
placed  in  an  area  of  image  2  that  contained  the  edge  device. 
A  second  ROI,  ROI2,  was  placed  in  an  area  of  image  2  that 
was  outside  of  the  edge  device.  The  detector  was  judged  to 
have  negligible  residual  signal  with  time  delay  r  if  it  passed 
the  following  test:45 


0  6  10  15  20  25  30 


(a) 


Exposure  (pC/kg) 


(b)  Exposure  (pC/kg) 

Fig.  3.  Plot  of  mean  pixel  value  versus  exposure  for  two  Mo/Mo  beams  and 
two  W/Rh  beams  over  the  (a)  entire  measured  exposure  range  and  (b)  the 
lower  exposure  range.  While  the  detector  exhibits  good  linearity  over  the 
entire  range,  divergences  from  linearity  occur  in  the  lower  exposure  range. 
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FlO.  4.  Plot  of  detector  MTF  along  (a)  CN  and  (b)  LR  axes  for  two  Mo/Mo 
beams  and  two  W/Rh  beams.  The  MTFs  for  the  four  beams  are  very  similar 
for  the  CN  axis,  but  differ  along  the  LR  axis.  The  pixel  aperture  limit  and 
theoretical  MTF  (Ref.  18)  are  included  for  reference. 


Table  III.  Summary  of  the  detector  MTF  properties  along  CN  and  LR  axes 
for  (top)  Mo/Mo  beams  and  (bottom)  W/Rh  beams.  Shown  are  frequencies 
at  specific  MTFs  and  the  MTF  at  specific  frequencies.  The  MTF  for  the 
Mo/Mo  beams  differed  between  the  CN  and  LR  axes,  but  was  similar  for 
the  W/Rh  beams. 


RQA-M2 
(CN  Axis) 

RQA-M2 
(LR  Axis) 

RQA-M4 
(CN  Axis) 

RQA-M4 
(LR  Axis) 

0.2  MTF 

11.1  mm-1 

12.7  mm-1 

1 1 .2  mm"1 

12.5  mm"1 

0.1  MTF 

12.8  mm’1 

14.8  mm"1 

12.8  mm"1 

14.5  mm"1 

0.5  mm-1 

0.98 

0.99 

0.99 

1.0 

2.5  mm-1 

0.86 

0.88 

0.89 

0.90 

5.0  mm-1 

0.65 

0.70 

0.67 

0.71 

MW2 

MW2 

MW4 

MW4 

(CN  Axis) 

(LR  Axis) 

(CN  Axis) 

(LR  Axis) 

0.2  MTF 

11.1  mm-1 

11.4  mm"1 

1 1.2  mm"1 

1 1.5  mm"1 

0.1  MTF 

12.9  mm"1 

13.2  mm"1 

12.9  mm"1 

13.3  mm"1 

0.5  mm-1 

0.98 

0.99 

0.99 

0.99 

2.5  mm-1 

0.86 

0.87 

0.88 

0.89 

5.0  mm-1 

0.65 

0.67 

0.67 

0.68 

ining  the  lower  exposure  range,  such  as  that  shown  in  Fig. 
3(b),  where  some  deviations  from  linearity  were  seen. 

To  verify  image  repeatability  over  time,  an  ensemble  of 
images  was  acquired  at  identical  mAs.  The  mean  signal  for 
each  image  was  computed  as  the  average  pixel  value  over  a 
region  of  interest.  These  images  showed  very  similar  signal 
levels  over  time,  as  the  mean  signal  varied  by  0.009%  over 
the  entire  ensemble  of  images.  In  contrast,  the  spatial  devia¬ 
tion,  which  described  how  the  pixel  values  varied  across 
each  image,  reached  3.4%  for  the  lowest  exposure  images. 

The  resolution  properties  of  the  detector,  as  represented 
by  the  MTF,  are  shown  in  Fig.  4  and  summarized  in  Table 
III.  Figure  4(a)  illustrates  the  MTF  along  the  CN  axis,  while 
Fig.  4(b)  displays  the  MTF  along  the  LR  axis.  While  the 
MTFs  for  the  tungsten  and  molybdenum  techniques  over¬ 
lapped  considerably  for  the  CN  direction,  the  molybdenum 


|(^|-^|)-(Cr3-^3)l 

\  +  % 

2 


0.005, 


(4) 


where  and  r)t  represented  the  mean  of  ROI,  and  ROI2  at 
time  t,  respectively.  The  IEC  chose  the  threshold  of  0.005  as 
the  maximum  allowable  level  of  residual  signal. 


III.  RESULTS 

Figure  3  illustrates  the  relationship  between  pixel  value 
and  exposure  for  the  detector.  In  general,  the  system  showed 
a  very  linear  response  with  correlation  coefficients  for  linear 
regression  fits  greater  than  0.999.  One  interesting  trend  was 
that  the  detector  was  slightly  more  sensitive  to  the  W/Rh 
beam  qualities  than  the  Mo/Mo  beams,  as  the  tungsten 
curves  resulted  in  higher  slopes  and  higher  pixel  values  for 
equivalent  exposures.  Another  trend  was  revealed  by  exam- 


-6  -4  -2  0  2  4  6 

Frequency  (mm'1) 


Fig.  5.  Two-dimensional  NNPS  for  RQA-M2  beam  quality  at  1.58  /rC/kg 
exposure.  The  image  is  shown  in  a  logarithmic  scale.  Nonstochastic  noise  is 
observed  in  a  frequency  band  along  the  CN  frequency  axis. 
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(b)  Frequency  (mm'1)  (d)  Frequency  (mm'1) 

Fig.  6.  Radial  traces  of  NNPS  multiplied  by  exposure  for  (a)  RQA-M2,  (b)  RQA-M4,  (c)  MW2,  and  (d)  MW4  beam  qualities. 


and  tungsten  curves  diverged  for  the  LR  axis.  For  reference, 
Fig.  4  also  displays  the  pixel  aperture  function  and  the  the¬ 
oretical  limit  calculated  by  Yorker  et  a/.18 

Figure  5  shows  an  example  of  a  two-dimensional  NNPS 
displayed  in  a  logarithmic  scale  (RQA-M2  technique, 
1.58  /zC/kg).  The  figure  demonstrates  nonstochastic  noise  in 
theCN  direction  along  a  band  of  0.1 12  mm-1  in  width.  Simi¬ 
lar  behavior  was  observed  for  other  exposures  and  beam 
qualities.  Figure  6  illustrates  the  radial  NNPS  multiplied  by 
exposure.  As  discussed  in  Sec.  II E,  the  product  of  NNPS 
and  exposure  should  remain  constant  for  strictly  quantum- 
limited  detectors,  however,  the  results  showed  notable  expo¬ 
sure  dependencies.  For  lower  exposures,  the  magnitude  of 
this  metric  decreased  to  some  minimum  value  as  one  in¬ 
creased  exposure.  For  several  techniques,  the  magnitude  of 
the  metric  increased  at  higher  exposures. 

Figure  7  illustrates  the  NNPS  calculated  through  the  back¬ 
ground  subtraction  method.  The  background  subtraction 
method  noticeably  reduced  the  low-frequency  noise.  In  addi¬ 
tion,  the  overall  magnitude  of  the  NNPS  decreased.  One 
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Fig.  7.  Radial  traces  of  NNPS  multiplied  by  exposure  obtained  with  and 
without  background  subtraction  method.  The  NNPS  was  obtained  using 
RQA-M2  technique  at  12.6  /zC/kg.  The  background  subtraction  routine  re¬ 
duced  the  low-frequency  noise  and  lowered  overall  noise. 
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Frequency  (mm'1) 
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Fig.  8.  DQE  averaged  over  CN  and  LR  axes  for  (a)  RQA-M2,  (b)  RQA-M4,  (c)  MW2,  and  (d)  MW4  beam  qualities. 


Table  IV.  Detector  DQE  along  the  (top)  CN  axis  and  (bottom)  LR  axes.  As  the  low-frequency  noise  caused 
peaks  in  the  DQE,  the  table  reports  the  maximum  DQE  value  and  the  frequency  at  which  this  maximum  occurs. 


RQA-M2 
(1.60  /tC/kg) 

RQA-M4 
(2.79  /uC/kg) 

MW2 

(1.90  jtC/kg) 

MW4 

(1.94  /uC/kg) 

Background  subtracted 
RQA-M2 
(12.6  fiC/kg) 

0.15  mm-1 

46% 

59% 

46% 

50% 

73% 

2.5  mm-1 

49% 

55% 

61% 

66% 

64% 

5.0  mm-1 

31% 

36% 

41% 

44% 

44% 

Peak 

55% 

63% 

66% 

77% 

73% 

1.25  mm-1 

0.85  mm-1 

1.55  mm'1 

0.85  mm-1 

0.15  mm'1 

Background  subtracted 

RQA-M2 

RQA-M4 

MW2 

MW4 

RQA-M2 

(1.60  /iC/kg) 

(2.79  /rC/kg) 

(1.90  /rC/kg) 

(1.94  /liC/kg) 

(12.6  /tC/kg) 

0.15  mm-1 

47% 

59% 

47% 

52% 

ITU 

2.5  mm”1 

50% 

57% 

61% 

70% 

68% 

5.0  mm-1 

34% 

41% 

41% 

47% 

49% 

Peak 

53% 

63% 

63% 

77% 

ITU 

1.45  mm"1 

0.95  mm"1 

1.85  mm-1 

0.85  mm 

0.15  mm'1 
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Fig.  9.  DQE  calculated  using  background  subtraction  method  averaged  over 
the  CN  and  LR  axes.  The  DQE  was  computed  for  RQA-M2  technique  at 
12.6  /rC/kg.  The  background  subtraction  routine  reduced  the  low  frequency 
peaking.  The  plot  also  shows  a  theoretical  estimation  of  the  DQE  (Ref.  46) 
for  a  similar  detector  (200  pm  selenium  layer,  85  pm  pixel  size). 

should  note  that  the  noise  was  corrected  for  the  change  in 
noise  variance  due  to  the  subtraction  technique,  as  discussed 
in  Sec.  II  E. 

Figure  8  shows  the  DQE  measured  for  all  techniques  in 
the  axial  (the  average  of  CN  and  LR  axes)  direction.  The 
DQE  curves  showed  low  frequency  peaking,  in  that  the  DQE 
exhibited  a  sharp  increase  at  lower  frequencies.  The  strong 
low-frequency  component  of  the  NNPS  led  to  this  unusual 
behavior.  Moreover,  the  DQE  increased  with  exposure  for 
lower  exposure  values,  reached  a  peak  value,  and  then  de¬ 
creased  for  higher  exposures.  This  was  also  expected  from 
the  behavior  of  the  NNPS.  The  DQE  is  summarized  in  Table 
IV  for  all  four  techniques  along  both  CN  and  LR  axes. 

To  separate  the  fixed  pattern  noise  from  quantum  noise 
effects,  the  DQE  was  calculated  with  the  background  sub¬ 
traction  method.  Figure  9  illustrates  the  DQE  computed  with 
this  method  in  the  axial  direction.  By  eliminating  the  fixed 
pattern  noise,  the  low-frequency  peaking  in  the  DQE  was 
removed  and  the  overall  efficiency  increased.  This  figure 
also  includes  the  theoretical  DQE  calculated  for  a  similar 
detector  for  reference  (200  pm  selenium  layer,  85  pm  pixel 
size).46 

The  results  from  lag  measurements  are  summarized  in 
Table  V.  In  general,  the  image  lag  for  the  detector  passed  the 
test  established  by  IEC  62220-1.  An  interesting  phenomenon 
occurred  for  the  fourth  test  (75  pGy  exposure,  5  min  decay 
time).  The  residual  signal  level  was  unacceptably  high  for 
this  test,  even  though  a  similar  test  (75  pGy  exposure,  3  min 
decay  time)  produced  acceptable  levels  of  residual  signal. 

IV.  DISCUSSION 

Digital  mammography  has  begun  to  replace  screen-film 
systems  in  some  clinical  settings.  The  motivation  for  this 
change  includes  several  logistical  considerations,  such  as 


Table  V.  Lag  properties  of  the  detector.  The  lag  tests  were  executed  in  the 
order  shown  in  the  table,  with  the  top  three  rows  measuring  signal  retention 
after  a  3  min  decay  time,  then  a  gap  of  1 0  mins,  and  the  bottom  three  rows 
measured  signal  retention  after  a  5  min  delay  time.  The  metric  corresponded 
to  IEC  62220-1  with  values  less  than  0.005  acceptable  under  the  IEC  guide¬ 
lines. 


Test  Number 

Exposure 

(pGy) 

Decay  time 
(min) 

Metric 

Acceptable  residual 
signal? 

1 

75 

3 

0.002 

Yes 

2 

150 

3 

0.0048 

Yes 

3 

200 

3 

0.0044 

Yes 

Ten  minute  wait 

4 

75 

5 

0.047 

No 

5 

150 

5 

0.013 

No 

6 

200 

5 

0.0022 

Yes 

convenient  archiving  and  display,  and  potential  image  quality 
advantages.  The  two  flat-panel  technologies  currently  of¬ 
fered,  direct  and  indirect,  vary  markedly  in  terms  of  their 
image  quality  characteristics.  Direct  detectors  tend  to  enjoy 
higher  resolution  than  indirect  detectors.  However,  they  are 
often  less  efficient  than  their  indirect  counterparts.9,10  In  this 
study,  we  evaluated  all  physical  properties  of  a  particular 
direct  flat-panel  detector,  including  resolution,  noise,  and  ef¬ 
ficiency,  to  enable  a  thorough  comparison  between  that  de¬ 
tector  and  others. 

Several  other  investigators  have  examined  the  physical 
characteristics  of  flat-panel  mammographic  imagers.  As 
such,  the  results  from  this  system  characterization  must  be 
reported  in  the  context  of  the  performance  of  other  systems. 
When  considering  previous  measurements,  one  should  note 
any  differences  in  beam  energies  and  filtrations.  Most  prior 
studies  utilized  molybdenum  anodes  with  molybdenum  fil¬ 
tration  at  28  kVp,  but  often  utilized  a  breast  equivalent  phan¬ 
tom  for  further  filtration.18,46  49  While  this  should  still  allow 
for  reasonable  comparisons  between  MTFs,  these  differences 
would  make  comparisons  of  DQE  curves  more  challenging. 

Compared  to  previous  measurements  of  indirect  flat-panel 
imagers,47’48  the  current  system  exhibited  a  higher  MTF.  At 
low  frequencies,  our  MTF  was  similar  to  other  direct  flat- 
panel  imagers,  but  our  MTF  was  higher  at  higher 
frequencies.46'49  As  the  pixel  size  served  as  the  primary  lim¬ 
iter  of  the  resolution  of  a  direct  detector,  with  some  blurring 
effects  from  backscatter  and  reabsorption  of  K  x-rays,  the 
similarity  between  direct  detectors  was  reasonable.  At  simi¬ 
lar  exposures,  the  DQE  of  the  system  was  generally  higher 
than  that  of  indirect  flat-panel  imagers,  although  the  low- 
frequency  peaking  complicated  this  comparison.47  In  com¬ 
parison  to  the  work  by  Jee,  the  High  Light  (HL)  output  con¬ 
figuration  produced  a  higher  DQE  but  the  High  Resolution 
(HR)  configuration  appeared  to  produce  a  lower  DQE  than 
our  current  system.48  The  direct  detector  evaluated  by  Zhao 
produced  a  generally  higher  DQE,  with  constant  behavior  for 
different  exposures,  although  there  were  significant  differ¬ 
ences  between  the  axes.46  One  interesting  result  was  that 
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Fig.  10.  Example  image  of  uniform  beam  of  radiation  (a)  before  and  (b)  after  the  application  of  a  secondary  gain  correction  from  the  average  of  10  images. 
The  larger  images  (physical  size:  23.3  cm  X  28.7  cm)  show  the  differences  in  large-scale  gain  nonuniformities.  Zoomed  portion  of  the  images  (2.1  cm 
X2.1  cm)  highlighting  pixel  artifacts  (c)  before  and  (d)  after  gain  calibration.  The  gain  calibration  largely  removes  the  pixel  artifacts  from  the  individual 
images. 
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Fig.  11.  Radial  traces  of  NNPS  multiplied  by  exposure  obtained  with  and 
without  background  subtraction  method  for  a  second  prototype  detector.  The 
NNPS  was  obtained  using  RQA-M2  technique  at  7.67  fiC/ kg.  As  with  the 
prototype  system,  this  detector  unit  also  exhibits  significant  stochastic  noise. 

some  previous  studies  on  selenium  detector  systems46,49  have 
also  observed  some  low-frequency  peaking  in  the  DQE,  al¬ 
though  it  was  less  pronounced  than  that  observed  in  this 
study. 

Yorker  et  al.  have  published  the  MTF  and  DQE  measure¬ 
ments  of  a  similar  mammographic  detector  with  an  identical 
pixel  size.18  That  study  examined  the  MTF  and  DQE  for  one 
radiographic  technique  using  a  molybdenum  anode  operated 
at  28  kVp  with  molybdenum  filtration.  The  reported  MTF  for 
this  technique  was  very  similar  to  our  measured  MTF  ac¬ 
quired  at  the  RQA-M2  technique.  At  similar  exposures,  our 
DQE  acquired  at  RQA-M2  also  appeared  comparable  to  that 
of  Yorker  et  al.  However,  given  the  fact  that  Yorker  et  al. 
used  a  4.2  cm  breast  phantom  filter  in  the  beam,  quantitative 
comparisons  are  not  straightforward. 

Several  researchers  have  explored  the  theoretical  proper¬ 
ties  of  selenium-based  flat-panel  imagers  operated  at  mam¬ 
mographic  energies.  The  properties  of  this  system  compared 
favorably  with  these  theoretical  calculations.  As  shown  in 
Figs.  4(a)  and  4(b),  the  MTF  of  this  system  remained  close  to 
the  theoretical  limit,  as  calculated  by  Yorker  et  a/.18  When 
using  the  background  subtraction  method  to  remove  fixed 
pattern  noise,  the  DQE  of  the  system  appeared  similar  to  its 
theoretical  value  for  a  similar  detector.46  The  difference  be¬ 
tween  the  theoretical  and  experimental  values  were  likely 
due  to  the  assumptions  behind  the  theoretical  calculation, 
which  assumed  a  200  pm  selenium  layer  and  85  pm  pixel 
size.  These  theoretical  calculations  should  underestimate  the 
actual  detector  efficiency,  as  a  larger  selenium  layer  will 
more  efficiently  capture  x-ray  photons  and  a  smaller  pixel 
size  should  boost  the  higher  frequency  portions  of  the  DQE. 
Notwithstanding,  the  experimental  results  for  the  MTF  and 
DQE  largely  agreed  with  their  theoretical  values. 

This  prototype  detector  had  very  favorable  resolution 
properties,  as  shown  by  its  MTF.  There  was  an  asymmetry  in 
the  MTF,  however,  as  the  tungsten  and  molybdenum  curves 


overlapped  for  the  CN  axis  but  diverged  for  the  LR  axis.  The 
difference  might  be  attributed  to  the  differences  between  the 
focal  spots  for  the  two  anodes,  in  terms  of  both  shape  and 
location.  The  impact  of  focal  spot  blur  should  be  minimal  as 
the  edge  was  placed  directly  on  the  detector.  Further  work 
remains  needed  to  evaluate  the  focal  spot  properties  for  both 
anode  varieties  and  to  determine  whether  this  phenomenon 
occurs  with  different  tubes. 

The  prototype  detector  did  show  significant  structured 
noise  contributions.  This  could  be  decomposed  into  two  fac¬ 
tors:  (a)  low  frequency  trending  over  the  image  and  (b)  pixel 
artifacts.  The  trending  was  expressed  as  a  strong  low- 
frequency  component  of  the  NNPS.  In  contrast,  pixel  arti¬ 
facts  were  similar  to  delta  functions  and  elevated  all  frequen¬ 
cies  of  the  NNPS.  A  background  subtraction  method 
eliminated  both  of  these  factors,  so  one  was  unable  to  deter¬ 
mine  the  relative  magnitude  of  either  individually.  Therefore, 
when  one  compared  the  NNPS  calculated  using  the 
background-subtraction  technique  to  that  calculated  using 
standard  techniques,  one  noticed  a  decrease  in  the  low- 
frequency  noise  contribution  as  well  as  an  overall  decrease  in 
the  magnitude  of  the  NNPS.  This  was  reflected  in  the  DQE 
as  well.  When  the  DQE  was  calculated  using  background- 
subtraction  techniques,  the  sharp  low-frequency  drop  was 
eliminated  and  the  overall  curve  was  shifted  upwards  be¬ 
cause  of  the  decrease  in  noise. 

Several  of  these  noise  concerns  could  be  mitigated  by 
additional  gain  calibration  after  the  gain  calibration  per¬ 
formed  by  the  system.  To  examine  the  benefits  of  further 
gain  calibration,  a  gain  map  was  created  by  averaging  ten 
uniform  images  together.  This  gain  map  was  then  applied  to 
a  subsequently  acquired  image.  The  effect  of  the  gain  cali¬ 
bration  is  shown  in  Fig.  10  and  displayed  with  identical  win¬ 
dow  and  level  settings.  The  prominent  trending  was  greatly 
diminished  and  many  of  the  pixel  artifacts  were  eliminated. 
To  assess  whether  the  problem  observed  was  unique  to  the 
prototype  detector  tested,  a  follow-up  experiment  was  con¬ 
ducted  on  a  more  recent  prototype  device  to  learn  whether  it 
exhibited  noise  properties  similar  to  the  earlier  prototype. 
This  experiment  compared  the  NNPS  calculated  with  and 
without  the  background  subtraction  methodology,  as  shown 
in  Fig.  11.  The  background  subtraction  proved  to  similarly 
remove  significant  nonuniformities,  which  indicated  that  the 
images  after  system  gain  calibration  retained  substantial 
structured  noise  in  the  second  prototype  as  well. 

The  detector  was  evaluated  for  four  different  beam  quali¬ 
ties.  Two  beams  used  a  molybdenum  anode  with  molybde¬ 
num  filtration  and  two  used  a  tungsten  anode  with  rhodium 
filtration.  The  tungsten  beam  qualities  were  developed  spe¬ 
cifically  for  this  study  and  inspired  by  IEC  standards.  The 
detector  appeared  to  be  slightly  more  sensitive  to  the  tung¬ 
sten  beams,  as  shown  in  the  exposure-pixel  value  relation¬ 
ship.  Moreover,  the  DQEs  for  the  tungsten  beams  were 
higher  than  those  for  the  molybdenum  beams,  although  this 
was  obscured  by  the  peaking  in  the  DQE  curves.  This  sug¬ 
gests  that  tungsten  beams  might  produce  higher  quality  im¬ 
ages  with  digital  detectors  than  the  traditional  molybdenum 
beams. 
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Using  the  parameters  established  by  IEC  62220-1,  the  im¬ 
age  lag  appeared  within  reasonable  parameters.  Nevertheless, 
high  exposures  led  to  unusual  behavior  in  signal  retention, 
affecting  other  exposures  even  after  a  significant  length  of 
time.  Ten  minutes  before  the  fourth  lag  test  (75  pGy  expo¬ 
sure,  5  min  decay  time),  a  200  iiGy  lag  test  was  conducted.  It 
appeared  that  this  high  exposure  still  affected  the  detector 
after  10  mins,  as  a  75  pGy  exposure  should  not  have  higher 
residual  signal  after  a  5  min  decay  time  than  it  would  after  a 
3  min  decay  time.  The  mechanisms  for  this  behavior  are 
unknown  and  suggest  additional  investigation  into  signal  re¬ 
tention  properties  of  selenium. 

V.  CONCLUSIONS 

This  study  reported  an  assessment  of  image  quality  for  a 
prototype  mammographic  imager  based  on  a  direct- 
detection,  flat-panel  array  employing  an  amorphous  selenium 
converter.  The  results  indicated  that  the  detector  had  strong 
potential  for  capturing  high-frequency  information,  as  exhib¬ 
ited  by  its  high  MTF.  In  addition,  the  DQE  of  the  detector 
approached  the  high  value  of  75%-80%.  Yet,  suboptimal 
calibration  affected  the  DQE  performance  of  the  system,  un¬ 
derscoring  the  importance  of  careful  gain  and  dead  pixel  cor¬ 
rections  in  reducing  detector  nonuniformities.  Finally,  this 
study  introduced  two  new  radiographic  techniques  utilizing 
tungsten  anodes  for  the  assessment  of  mammographic  sys¬ 
tems,  which  will  facilitate  the  future  comparisons  of  detector 
characteristics  operated  with  tungsten  anodes. 
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ABSTRACT 

For  diagnosis  of  breast  cancer  by  mammography,  the  mammograms  must  be  viewed  by  a  radiologist.  The  purpose 
of  this  study  was  to  determine  the  effect  of  display  resolution  on  the  specific  clinical  task  of  detection  of  breast 
lesions  by  a  human  observer.  Using  simulation  techniques,  this  study  proceeded  through  four  stages.  First,  we 
inserted  simulated  masses  and  calcifications  into  raw  digital  mammograms.  The  resulting  images  were  processed 
according  to  standard  image  processing  techniques  and  appropriately  windowed  and  leveled.  The  processed  images 
were  blurred  according  to  MTFs  measured  from  a  clinical  Cathode  Ray  Tube  display.  JNDMetrix,  a  Visual 
Discrimination  Model,  examined  the  images  to  estimate  human  detection.  The  model  results  suggested  that 
detection  of  masses  and  calcifications  decreased  under  standard  CRT  resolution.  Future  work  will  confirm  these 
results  with  human  observer  studies.  (This  work  was  supported  by  grants  NIH  R21-CA95308  and  USAMRMC 
W  8 1 XWH-04- 1-0323.) 

Keywords:  Image  Quality,  Mammography,  Simulation,  Task-Based  Assessment 


1.  INTRODUCTION 

After  a  digital  mammogram  has  been  acquired,  a  human  observer  must  view  the  data  in  order  to  detect  or  diagnose 
disease.  The  display  device,  therefore,  assumes  a  crucial  role  in  the  imaging  chain.  While  several  researchers  have 
given  significant  attention  to  the  quality  of  image  acquisition,1'9  fewer  investigators  have  measured  the  impact  of 
display  devices.10 13  To  understand  this  impact,  studies  must  evaluate  the  physical  properties  of  these  devices. 
However,  while  physical  characterization  remains  important,  display  quality  must  ultimately  be  described  in  terms 
of  the  clinical  task  in  question.14  16  This  study  considered  this  type  of  question,  examining  the  impact  of  display 
resolution  on  the  detection  of  mammographic  lesions. 

A  Cathode  Ray  Tube  (CRT)  display  serves  as  a  common  mammographic  display  device.17  As  a  CRT  ages,  its 
resolution  becomes  progressively  more  degraded,  leading  to  lower  display  quality  over  time.18  The  purpose  of  this 
study  was  to  consider  how  this  degradation  in  resolution  impacted  the  clinical  utility  of  a  CRT  display,  specifically 
the  detection  of  breast  masses  and  calcifications. 

2.  METHODS  AND  MATERIALS 

In  this  study,  first  simulated  masses  and  calcifications  were  inserted  into  digital  mammograms.  We  applied  basic 
image  processing  techniques  to  these  images  and  adjusted  the  window  and  level  appropriately.  Next,  we  blurred  the 
images  according  to  three  different  resolution  settings  measured  from  a  CRT  display.  Finally,  a  model  observer 
viewed  each  of  these  images  to  estimate  the  detection  probabilities  under  each  blur  setting.  The  following  describes 
the  details  of  these  steps. 

2.1  Acquisition  of  Digital  Mammographic  Backgrounds 

Digital  mammographic  images  were  acquired  on  a  clinical  flat-panel  cesium  iodide-based  digital  mammography 
system  (Senographe  2000D,  GE  Medical  Systems,  Milwaukee,  WI).  Previous  studies  have  characterized  the 


physical  characteristics  of  this  digital  mammography  system.19,  20  Images  used  in  this  study  were  normal 
craniocaudul  view  mammograms  acquired  with  a  molybdenum  anode  with  molybdenum  or  rhodium  filtration.  The 
beam  energies  for  the  images  ranged  from  25  to  30  kVp  and  compressed  breast  thicknesses  extended  from  2.7  cm  to 
7.3  cm  with  varying  glandular  and  adipose  tissue  composition. 

2.2  Lesion  Simulation 

Simulated  breast  lesions  were  placed  in  the  center  of  mammographic  images  using  an  established  procedure  for 
simulating  masses  and  calcifications  with  attributes  similar  to  those  of  real  mammographic  lesions.21, 2  Breast  mass 
simulation  proceeded  through  three  stages,  as  illustrated  in  Figure  1 .  The  first  stage  sets  each  pixel  of  an  array  to  its 
equivalent  major  axis  value, 

b  =  I(y  -  yo)Cos  [a ] -  (x -xg)Sin  [a]]2  +  [(* -jr0)Cas  [a]  +  (y -  yQ)Sin  [a ]f  W 


where  (xq,  y0)  represent  the  center  of  the  mass,  a  determines  the  angular  orientation  of  the  mass,  and  c  corresponds 
to  the  ratio  of  the  minor  axis  length  to  the  major  axis  length.  The  second  stage  introduced  non-uniformities  in  the 
mass  border  by  multiplying  the  elliptical  rings  with  a  border  deviation  profile  with  a  given  variance  and  power 
spectrum.  The  final  stage  converted  the  equivalent  major  axis  values  to  detector  gray  level  values  through  the 
elliptical  trace  function. 

The  calcification  procedure  similarly  required  three  stages.  The  first  stage  established  the  distribution  of 
calcifications,  using  either  a  clustered  or  linear  distribution.  The  second  stage  created  individual  calcification  at 
each  point  specified  by  the  calcification  distribution  through  a  series  of  morphological  thickening  and  erosion 
operations.  This  resulted  in  a  binary  mask  of  the  calcifications.  The  final  stage  added  the  binary  mask  to  a 
background  image  with  the  appropriate  contrast. 

The  spatial  parameters  for  the  simulation  routine  were  determined  from  screen-film  mammographic  data  obtained 
through  the  Digital  Database  of  Screening  Mammography.23  These  parameters  remained  applicable  to  digital 
mammographic  backgrounds.  However,  the  lesion  contrast  must  be  separately  calculated  for  the  digital  case  as  the 
contrast  in  screen-film  images  were  impacted  by  varying  H&D  characteristics.  To  determine  the  appropriate 
contrast  for  the  simulated  lesions,  the  xSpect  x-ray  simulation  program24  calculated  the  unit  contrast  for  both  masses 
and  calcifications  embedded  in  a  50%  glandular/50%  adipose  breast  imaged  with  a  cesium-iodide  detector.  The 
contrasts  were  calculated  for  a  molybdenum  anode  with  molybdenum  or  rhodium  filtration  for  each  kVp  and  breast 
thickness.  Contrast  reduction  by  scattered  radiation  was  also  accounted  for  using  previously  published 
measurements.25  The  lesions  were  then  inserted  in  mammographic  backgrounds  with  the  appropriate  contrast  and 
spatial  features  for  the  given  anode,  filtration,  kVp,  and  breast  thickness. 
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Fig.  1  Schematic  of  mass 
simulation  procedure.  The 
three  images  illustrate  the  three 
steps  in  this  system. 
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Fig.  2  Illustration  of 
calcification  simulation 
routine.  This  procedure 
considers  each  point  in 
the  distribution  and 
creates  a  unique 
calcification  on  each 
point. 


2.3  Image  Processing 

Most  digital  mammography  systems  employ  post-processing  algorithms  to  improve  image  display.  A  common 
technique  separates  the  images  into  multiple  frequency  bands  to  improve  contrast  for  specific  frequencies.  This 
study  utilized  a  basic  image  processing  algorithm  that  enhanced  two  frequency  bands  in  the  image.26  The  first  stage 
augmented  the  higher  frequency  content  of  the  image,  while  the  second  stage  strengthened  the  content  variations. 
The  parameters  for  each  stage  were  determined  by  visual  analysis  of  the  images.  The  first  stage  accentuated  the 
sharp  detail  in  the  image  through  an  unsharp  masking  procedure  as, 

lvs  =/  +  5F(c)-(/-E@/)  (2) 

where  /  represented  the  input  image,  L,  the  Gaussian  kernel,  had  a  standard  deviation  of  0.45  mm  and  width  of  2.8 
mm,  and  SF(c),  the  sharpness  factor,  controlled  the  level  of  enhancement.  To  boost  low  contrast  objects,  a  non¬ 
linear  function  was  utilized  for  SF(c)  as, 


SF(c)  = 


G  1- 


c0 

SF  o 


+  SF0  c<c0 


C  >  Cq 


(3) 


with  a  gain,  G,  of  1,  a  contrast  threshold,  c0,  of  70,  a  contrast,  c,  equal  to  the  absolute  difference  between  the  blurred 
image  and  original  image,  and  a  slope  parameter,  p,  of  3.  The  second  stage  enhanced  the  mid-frequency 
components  of  the  image,  as 

l  On,  =(Z  ®  I  us) +  CF(c )  •  (Ius  -£®IUS)  (4) 

where  £  represented  a  Gaussian  kernel  with  a  standard  deviation  of  4.4  mm  and  CF(c)  controlled  the  level  of 
contrast  enhancement.  The  function  CF(c)  had  the  same  functional  form  as  SF(c),  but  utilized  a  gain,  G,  of  1.3. 


Once  the  images  were  processed,  observers  window  and  level  an  image  in  order  to  produce  an  acceptable  image 
appearance.  To  determine  the  window  and  level  parameters  for  each  mammogram,  an  experienced  mammographer 
windowed  and  leveled  each  mammogram  individually.  A  sigmoid  transformation  was  fit  to  each  window  and  level 
function,  to  provide  a  smooth  transition  at  the  extremes  of  the  display  pixel  values.  This  transformation  was 
represented  as 


hind  a[\  l+e-(Log[r0ut]-S  )/<r)  (5) 

where  I0ui  represented  the  processed  mammogram,  8 equaled  the  center  of  the  sigmoid  transition,  and  a  established 
the  slope  of  the  sigmoid  transition. 

2.4  Measurement  of  Display  Characteristics 

We  measured  the  resolution  properties  of  a  five  mega-pixel  Cathode  Ray  Tube  (CRT)  display  system  (Barco  MGD- 
521,  p45  phosphor)  with  a  10-bit  graphics  controller  (Barco  5MP1H).  A  charge-coupled  device  (CCD)  camera 
(XCD-SX900,  Sony  Corporation,  Tokyo,  Japan)  equipped  with  a  macro  lens  (Rodgen  1:4,  28mm,  Rodenstock, 
Munich,  Germany)  acquired  images  of  line  test  patterns  presented  on  the  display.  Two  images  from  the  recent 
TG18  test  pattern  set  (TG18-RV50  and  TG18-RH50)  supplied  a  vertical  and  horizontal  line,  respectively.27, 28  To 
remain  in  the  quasi-linear  range  of  the  display,  these  patterns  employed  subtle  lines,  with  12%  contrast  from  the 
background.  We  then  computed  the  MTF  from  these  line  patterns  using  established  methods.  Full  details  of  the 
measurement  methodology  has  been  reported  in  another  publication.29  We  measured  the  MTF  for  the  standard 
display  resolution  setting  and  two  degraded  resolution  settings  using  the  defocusing  feature  of  the  display.  These 
measured  MTFs  are  displayed  in  Figure  3. 


FlG.  3  Measured  MTF  for  a 
CRT  display  under  three  different 
resolution  settings. 


2.5  Simulation  of  Image  Display 

A  Resolution  Modification  routine,  the  details  of  which  are  disclosed  in  a  previous  publication,30  simulated  the  blur 
effects  of  the  CRT  display.  This  routine  altered  the  resolution  of  an  input  image  according  to  an  input  MTF  to 
produce  a  blurred  version  of  the  image.  To  accomplish  this,  the  input  mammogram  was  transformed  to  the 
frequency  domain  through  an  FFT.  The  frequency  content  of  the  image  was  then  filtered  by  the  display  MTF.  An 
inverse  FFT  transformed  this  modified  frequency  spectrum  back  to  the  spatial  domain.  This  blurring  was  performed 
for  each  display  resolution  setting  to  produce  multiple  versions  of  each  image. 


2.6  Observer  Model  Experiment 

A  5.12  cm  x  5.12  cm  region  of  interest  (ROI)  was  extracted  from  the  central  breast  area  for  analysis  by  a  visual 
discrimination  model  (VDM).  The  Samoff  JNDmetrix33  VDM  has  been  used  to  simulate  the  effects  of  display 
characteristics  and  image  processing  on  the  conspicuity  of  mammographic  lesions.11, 31, 32,34  For  this  study,  the 
VDM  compared  a  mammogram  containing  a  lesion  to  the  same  mammogram  without  the  lesion  and  computed  a 
just-noticeable  difference  (JND)  metric  for  the  discriminability  of  those  images  by  a  human  observer.  The  VDM 
first  convolved  the  input  images  by  an  approximation  of  the  point-spread  function  of  the  optics  of  the  eye.  The 
model  simulated  sampling  of  the  image  by  retinal  cones  by  performing  a  Gaussian  convolution  and  then  point 
sampling.  Next,  it  computed  the  local  contrast  from  the  raw  luminance  image.  The  model  applied  a  Laplacian 
pyramid  to  the  data  in  order  to  isolate  five  frequency  bands  from  the  data.  For  each  frequency  band,  the  data  was 
convolved  with  eight  pairs  of  spatially  oriented  filters.  The  sensitivities  and  other  parameters  for  these  filters  were 
determine  by  fitting  model  output  to  psychophysical  data  from  sine-grating  detection  and  discrimination 
experiments.  The  model  squared  each  pair  of  filter  output  images  and  summed  them  to  provide  a  phase-independent 
response.  Next,  the  transducer  stage  derived  the  energy  for  each  frequency  band,  normalizing  this  energy  by  the 
square  of  the  appropriate  grating  contrast  detection  threshold.  A  sigmoid  function  was  applied  to  each  frequency 
bands  to  account  for  the  visual  contrast  discrimination  function.  The  model  incorporated  the  foveal  sensitivity  by 
averaging  the  outputs  from  the  transducer  step  using  a  disk  kernel.  The  final  product  of  the  model  was  a  two- 
dimensional  map  of  JND  values,  where  each  pixel  indicated  the  discriminability  of  the  two  input  images. 

3.  RESULTS 

Figure  4  illustrates  the  results  when  the  VDM  discriminates  between  mammographic  images  with  simulated  benign 
masses  and  those  without  simulated  benign  masses.  The  perfect  resolution  refers  to  images  without  any  display 
blur,  while  the  other  three  resolution  settings  refer  to  the  measured  MTFs  in  Figure  3.  As  expected,  the  model  was 
better  able  to  detect  masses  without  any  display  blur.  The  difference  between  the  three  display  blur  settings 
remained  much  more  modest. 

Fig.  4  JND  Aggregate  Measure 
(JAM)  from  the  VDM  comparing 
mammographic  backgrounds  with 
and  without  simulated  benign 
masses  for  four  different  resolution 
settings.  The  benign  masses  had 
an  average  diameter  of  3  mm.  The 
error  bars  represent  the  95% 
Confidence  Interval. 
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Figure  5  illustrates  the  results  when  the  VDM  discriminates  between  mammographic  images  with  simulated  fine 
linear  branching  calcifications  and  those  without  simulated  calcifications.  The  nomenclature  in  Figure  5  remains 
consistent  with  Figure  4.  As  expected,  the  model  had  a  greater  ability  to  detect  these  calcifications  without  any 
display  blur.  However,  as  the  resolution  of  the  CRT  degrades,  the  detectability  of  calcifications  decreased 
significantly.  Similar  model  results  were  obtained  for  images  with  pleomorphic  calcifications. 
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Fig.  5  JND  Aggregate  Measure 
(JAM)  from  the  VDM  comparing 
mammographic  backgrounds  with 
and  without  simulated  fine  linear 
branching  calcifications  for  four 
different  resolution  settings. 


4.  DISCUSSION  AND  CONCLUSIONS 

This  study  evaluated  the  impact  of  display  blur  on  the  detection  of  mammographic  masses  and  calcifications.  As  an 
initial  step,  this  study  utilized  a  visual  discrimination  model  to  estimate  detection  by  a  human  observer.  These  initial 
results  suggested  that  detection  of  masses  and  calcifications  decreased  with  standard  CRT  resolution.  In  addition, 
the  model  implies  that  the  detection  of  calcifications,  but  not  masses,  declined  as  the  resolution  of  the  CRT  degraded 
over  time.  This  prediction  seems  reasonable  because  the  conspicuity  of  small,  fine  structures  in  calcifications  are 
more  likely  than  larger  objects,  such  as  masses,  to  be  affected  by  reductions  in  display  resolution.  The  next  phase  of 
modeling  will  use  VDM  output  to  predict  signal  detectability  within  the  framework  of  a  channelized  model 
observer.  Future  work  must  include  human  observer  performance  experiments  to  verify  these  estimates. 
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Abstract 

Lesion  simulation  provides  a  tool  when  quantifying  the  utility  of  an  imaging  system  in  a 
detection  task.  For  mammography,  the  important  detection  tasks  are  detecting  breast  masses  and 
calcifications.  In  this  study,  we  characterized  the  radiographic  appearance  of  both  masses  and 
calcifications  from  images  obtained  from  the  Digital  Database  of  Screening  Mammography 
(DDSM).  The  characterization  results  were  then  used  in  a  routine  capable  of  creating  simulated 
masses  and  calcifications.  To  verify  the  quality  of  this  simulation  routine,  an  observer 
performance  experiment  was  conducted  in  which  an  observer  was  asked  to  discriminate  between 
real  and  simulated  lesions.  The  results  were  then  analyzed  using  ROC  analysis.  The  preliminary 
results  showed  an  Az  of  0.59  for  benign  masses,  0.61  for  malignant  masses,  and  0.58  for 
malignant  calcifications.  More  observer  studies  are  underway  to  enhance  the  statistical  power  of 
these  results.  (This  work  was  supported  by  a  grant  from  the  NIH,  R21-CA95308  and 
USAMRMC  W8 1 XWH-04- 1  -0323 .) 


1.  Introduction 

A  number  of  new  full-field  digital  mammography  systems  with  varying  attributes  have  entered 
the  clinical  arena.  It  is  important,  therefore,  to  discover  which  systems  are  most  appropriate  for 
mammographic  imaging.  As  the  detection  of  breast  cancer  is  the  key  task  in  mammography,  a 
system  should  be  judged  in  how  well  it  aids  in  that  task.  Simulation  techniques  significantly 
facilitate  such  evaluations  for  a  variety  of  detectors,  breast  densities,  and  lesion  types. 


One  hurdle  faced  by  mammography  simulation  is  the  lack  of  breast  lesion  models.  For  masses, 
previous  work  has  used  gaussian  profiles,  disks,  and  simulated  lung  nodule  profiles. 14  For 
calcifications,  the  most  common  model  has  been  to  utilize  masks  extracted  from  real 
calcifications.5,6  This  study  adopted  a  different  approach.  First,  we  characterized  the 
radiographic  appearance  of  breast  masses  and  calcifications  from  real  mammograms.  Then,  we 
created  simulated  breast  masses  and  calcifications  emulating  those  characteristics.  Our  mass 
model  was  previously  validated  through  a  preliminary  observer  performance  experiment.  This 
paper  extends  that  work  to  microcalcifications. 

2.  Lesion  Characterization 

2. 1  Breast  Mass  Characterization  Procedure 

Four  categories  of  breast  masses  were  chosen  for  characterization  using  the  BI-RADS  lexicon. 
Two  types  were  typically  benign,  oval  circumscribed  and  oval  obscured  masses,  and  two  were 
typically  malignant,  irregular  ill-defined  and  irregular  spiculated.  Sample  mammograms 
containing  these  lesions  were  extracted  from  the  University  of  South  Florida’s  Digital  Database 
for  Screening  Mammography  (DDSM).9  Characterization  was  performed  on  a  2.56  cm  x  2.56 
cm  region  of  interest  (ROI)  surrounding  the  mass.  All  ROIs  were  converted  to  optical  density 
values  using  the  characteristic  curve  of  the  scanner. 

The  behavior  of  the  masses  was  determined  through  a  large-scale  analysis  and  a  small-scale 
analysis.  The  large-scale  behavior  was  characterized  through  an  elliptical  trace,  which  examined 
the  changes  in  optical  density  through  concentric  elliptical  rings.  The  small-scale  behavior  was 
measured  through  a  deviation  profile  that  measured  how  the  border  of  the  lesion  varied  from  an 
ellipse.  These  are  shown  graphically  in  figure  1 . 


Fig.  1.  The  elliptical  trace,  left,  characterizes  the  large-scale  behavior  of  the  mass.  The  small- 
scale  behavior  is  shown  in  the  border  deviation  profile,  right. 

2.2  Breast  Mass  Characterization  Results 

Example  characterization  results  for  typically  benign  masses  are  shown  in  figure  2.  The 
elliptical  trace  showed  a  sharp  transition  from  the  mass  to  the  background,  which  was  expected 
for  a  circumscribed  border.  The  border  deviation  profile  showed  some  deviations  from  the 
perfect  ellipse,  but  the  magnitude  was  fairly  small.  This  was  in  contrast  to  the  results  for 
typically  malignant  masses,  an  example  of  which  is  shown  in  figure  3.  The  elliptical  trace  for 
these  masses  showed  a  very  slow  transition  from  the  mass  to  the  background.  The  border 
deviation  profile  illustrated  strong  deviations  from  the  perfect  ellipse.  This  was  expected  as  the 
borders  are  ill-defined  and  the  shape  was  irregular. 


Fig.  2.  Example  characterization  results  for  benign  masses.  The  elliptical  trace,  left,  shows  a 
strong  transition  from  mass  to  background  while  the  border  deviation  profile,  right,  shows  small 
deviations  from  the  perfect  ellipse. 


Equivalent  Major  Axis  (mm)  Angle  (Radians) 


Fig.  3.  Example  characterization  results  for  malignant  masses.  The  elliptical  trace,  left,  shows  a 
smooth  transition  to  background,  and  the  border  deviation  profile,  right,  shows  marked 
deviations  from  the  perfect  ellipse. 

2.3  Calcification  Characterization  Procedure 

Two  categories  of  calcifications  were  chosen  based  on  the  BI-RADS®  lexicon.8  The  two 
categories  were  fine  linear  branching  and  pleomorphic,  referring  to  typically  malignant  lesions. 
The  distribution  studied  for  fine  branching  calcifications  was  linear,  while  the  distribution 
studied  for  pleomorphic  was  clustered.  Similar  to  masses,  sample  mammograms  were  drawn 
from  the  DDSM  database. 

To  characterize  the  calcifications,  a  mask  of  the  distribution  was  drawn.  Measurements  were 
then  made  on  this  binary  mask.  Three  properties  were  measured  for  the  individual  calcifications: 
the  major  axis,  minor  axis,  and  the  average  contrast.  Furthermore,  the  distributions  for  each 
calcification  type  were  measured.  For  pleomorphic  calcifications,  the  major  axis  and  minor  axis 
of  the  cluster  were  measured.  For  fine  linear  branching  calcifications,  the  lengths  of  the  lines 
were  measured  along  with  the  angle  between  the  lines  of  calcifications. 


2.4  Calcification  Characterization  Results 


The  results  from  the  calcification  characterization  are  shown  in  table  1 .  The  individual 
calcifications  results  were  similar  for  both  pleomorphic  and  fine  linear  branching  categories. 
The  distribution  results  established  the  mean  shape  for  each  distribution. 


Table  1 .  Summary  of  calcification  characterization  results 


Calcifications: 

Pleomorphic 

Fine  Linear  Branching 

Major  Axis  (mm) 

0.47  ±0.11 

0.43  ±0.13 

Minor  Axis  (mm) 

0.29  ±0.057 

0.26  ±  0.045 

Contrast 

0.22  ±0.13 

0.34  ±0.16 

Distribution: 

Major  Axis  (mm) 

8.0  ±3.5 

n/a 

Minor  Axis  (mm) 

7.1  ±3.2 

n/a 

Line  Length  (mm) 

n/a 

6.2  ±  2.3 

Angle  (degrees) 

n/a 

50.8  ±  11.2 

3.  Lesion  Simulation 

3.1  Mass  Simulation 

The  mass  simulation  routine  began  with  an  array  where  each  pixel  was  set  equal  to  its  equivalent 
major  axis  value  (given  the  eccentricity  and  center  location).  The  border  deviation  effects  were 
then  applied  to  this  array.  Finally,  the  array  was  transformed  to  optical  density  using  the 
elliptical  trace  profile.  This  is  shown  graphically  in  figure  4.  Example  masses  are  shown 
imbedded  in  backgrounds  in  figure  5. 


Fig.  4.  Graphical  overview  of  mass  simulation  procedure.  The  image  on  left  shows  an  array 
with  pixel  values  equal  to  their  equivalent  major  axis  value.  The  border  deviations  are 
introduced  in  the  center  image.  Finally,  the  image  is  transformed  to  optical  density  through  the 
elliptical  trace  profile,  which  results  in  the  final  image  on  the  right. 


Fig.  5.  Example  simulated  masses.  The  image  on  the  left  is  a  simulated  benign  mass,  while  the 
image  on  the  right  is  a  simulated  malignant  mass  with  an  ill-defined  border. 

3.2  Calcification  Simulation 

The  measured  distribution  results  established  a  probability  distribution  for  the  individual 
calcification  centers.  For  the  pleomorphic  category,  the  centers  had  a  uniform  probability 
density  within  an  ellipse  with  a  given  major  axis  and  minor  axis  length.  For  the  fine,  linear 
branching  case,  the  centers  had  a  uniform  probability  distribution  along  lines  with  a  given  mean 
length  and  relative  angle  between  lines. 

Given  the  desired  number  of  individual  calcifications,  the  simulation  program  sampled  these 
distributions  to  determine  the  location  of  the  centers  of  the  individual  calcifications.  For  each 
individual  calcification,  a  line  was  drawn  through  this  center  at  a  random  angle.  The  length  of 
this  line  was  equal  to  the  major  axis  length  of  the  individual  calcifications.  A  morphological 
thickening  operation  was  then  applied,  followed  by  a  morphological  eroding.  These  produced 
the  shapes  of  the  individual  calcifications.  The  calcification  distribution  was  then  added  to  a 
normal  background  with  a  given  contrast.  Example  simulated  calcifications  were  shown  in 


figure  6. 


Fig.  6.  Example  simulated  calcifications.  The  left  image  shows  a  simulated  pleomorphic 
distribution,  while  the  center  and  right  image  show  simulated  fine,  linear  branching 
calcifications. 


4.  Observer  Performance  Experiment 

4.1  Observer  Protocol 

To  determine  the  quality  of  the  simulation  routines,  an  observer  performance  experiment  was 
conducted.  In  this  study,  an  experienced  mammographer  was  asked  to  rate  their  confidence  in 
whether  a  lesion  was  definitely  real  or  definitely  simulated.  The  simulation  routine  would  be 
effective  if  a  mammographer  was  unable  to  distinguish  the  difference  between  the  simulated  and 
real  lesions.  As  this  was  a  preliminary  experiment,  only  one  mammographer  was  used. 


4.2  Observer  Results 


The  histograms  for  the  observer  results  for  masses  are  shown  in  figure  7.  In  general,  the 
distributions  for  real  and  simulated  masses  overlap  considerably.  The  histogram  for 
calcifications  is  shown  in  figure  8.  Again,  the  histograms  for  real  and  simulated  lesions  overlap 


considerably. 


Realism  Rating  Value 


Realism  Rating  Value 

Fig.  7.  Histograms  of  the  rating  frequency  versus  rating  value  for  real  and  simulated  masses. 
The  results  for  typically  benign  masses  are  shown  on  the  top  while  the  typically  malignant  mass 
results  are  shown  on  the  bottom.  The  typically  malignant  masses  are  further  separated  by  border 
type  for  real  and  simulated  masses. 


Realism  Rating  Value 

Fig.  8.  Histogram  of  the  rating  frequency  versus  rating  value  for  calcifications. 


To  quantify  the  degree  of  overlap,  a  Receiver  Operating  Characteristic  Analysis  was 
performed.10  In  this  case,  an  Az  of  0.5  indicates  that  an  observer  was  near  chance  in 
discriminating  between  real  and  simulated  lesions.  This  analysis  is  summarized  in  table  2. 


Table  2.  Summary  of  ROC  Analysis  for  discrimination  between  real  and  simulated  lesions. 


Az 

a 

Benign  Masses 

0.59 

0.08 

Malignant  Masses 

0.61 

0.07 

Malignant  Calcifications 

0.58 

0.07 

5.  Conclusions 

The  characterization  procedure  undertaken  in  this  study  introduces  a  new  way  to  describe  breast 
lesions.  The  data  from  this  characterization  was  then  used  in  a  new  simulation  routine  that  is 
capable  of  simulating  breast  masses  and  calcifications.  Results  from  a  preliminary  observer 
performance  experiment  on  these  simulations  indicate  that  our  simulation  routine  produces  high 
quality  simulations  of  breast  masses  and  calcifications.  Further  work  is  needed  to  validate  the 
results  of  this  preliminary  observer  performance  experiment. 
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Abstract 

Digital  mammography  has  the  potential  to  improve  image  quality  for  mammographic  imaging. 
This  study  evaluated  a  selenium-based  direct  full-field  digital  mammographic  imager  (70  (im 
pixels)  using  a  molybdenum  anode  operated  at  28  kVp  with  inherent  filtration  of  30  pm 
molybdenum  and  an  additional  2  mm  of  aluminum  filtration.  To  capture  the  detector  resolution, 
we  measured  the  presampled  modulation  transfer  function  (MTF)  using  an  edge  method.  The 
noise,  summarized  through  the  Normalized  Noise  Power  Spectrum  (NNPS),  was  measured  by 
two-dimensional  Fourier  analysis  of  uniformly  exposed  radiographs.  The  detective  quantum 
efficiency  (DQE)  was  then  computed  from  the  measured  MTF,  NNPS,  and  ideal  signal-to-noise 
ratio.  For  the  Left-Right  axis,  the  MTF  reached  the  value  of  0.2  and  0.1  at  12.7  mm'1  and  14.8 
mm"1,  respectively.  The  DQE  attained  a  maximum  value  of  53%  at  1 .45  mm'1  for  the  Left-Right 
axis.  However,  the  DQE  showed  a  strong  dependence  on  exposure  and  frequency.  The  results 
indicated  that  this  detector  has  high  resolution,  but  it  may  be  valuable  to  remove  structured  noise 
through  improved  calibration  before  clinical  implementation.  (The  full  data  for  this  study  are 
published  as  R.S.  Saunders,  Jr,  E.  Samei,  J.L.  Jesneck,  and  J.Y.  Lo,  "Physical  characterization  of 
a  prototype  selenium-based  full  field  digital  mammography  detector,"  Med.  Phys.  32(2)  (2005). 


1. 


Introduction 


The  purpose  of  this  work  was  to  evaluate  the  physical  characteristics  of  a  selenium  full-field 
digital  mammography  (FFDM)  detector.  Three  different  metrics  of  system  performance  were 
evaluated:  the  Modulation  Transfer  Function  (MTF),  Normalized  Noise  Power  Spectrum 
(NNPS),  and  Detective  Quantum  Efficiency  (DQE).  As  previous  research  has  shown  that 
selenium  detectors  can  exhibit  image  lag  and  ghosting,1  this  research  also  examined  the  lag 
performance  of  the  detector. 

2.  Methods  and  Materials 

2.1  Detector  Description 

The  detector  investigated  in  this  study  was  a  selenium-based  flat-panel  detector  (Siemens 
Medical  Systems,  Erlangen,  Germany).  The  detector  was  based  on  a  250  pm  selenium 
photoconductive  layer  coupled  to  a  storage  capacitor  and  amorphous  selenium  switching 
transistor.2  The  active  detector  area  was  23.296  cm  x  28.672  cm  consisting  of  3328  x  4096 
square  pixels  with  a  70  pm  pixel  pitch.  Prior  to  evaluation,  the  antiscatter  grid  supplied  with  the 
system  was  removed  and  gain  and  dead  pixel  corrections  were  performed  according  to 
manufacturer  specifications.  For  most  measurements,  the  standard  detector  cover  was  kept  in 
place  and  the  compression  paddle  was  removed.  For  the  MTF  measurements,  the  detector  cover 
was  removed  so  that  the  edge  device  could  be  placed  in  contact  with  the  detector. 

2.2  Image  Acquisition 

The  selenium  detector  was  coupled  to  a  high  frequency  multiphase  x-ray  generator  (Mammomat 
Novation)  for  which  the  high-voltage  accuracy  was  certified  to  be  within  ±5%.  All  images  were 


acquired  with  a  large  focal  spot  of  0.3  mm,  nominal.  We  used  the  RQA-M2  technique,  3,  which 
employed  a  molybdenum  anode  operated  at  28  kVp,  30  pm  molybdenum  inherent  filtration,  and 
2  mm  aluminum  added  filtration.  The  image  data  were  acquired  in  a  raw  format  without  any 
image  post-processing  applied.  After  acquisition,  the  images  were  transferred  to  a  research 
computer  as  14-bit,  raw  data  for  analysis. 

For  all  image  acquisitions,  the  exposure  to  the  detector  was  measured  free  in  air  using  a 
calibrated  ionization  chamber  (1515  x-ray  monitor  with  10X5-6M  dedicated  mammography 
ionization  chamber,  Radcal  Corporation,  Monrovia,  CA).  The  chamber  was  placed  17  cm  above 
the  detector  to  minimize  contributions  from  backscatter.  The  exposure  incident  on  the  detector 
at  65  cm  source  to  image  distance  (SID)  was  estimated  from  the  measured  exposure  using  the 
inverse-square  law. 

2.3  Linearity 

The  linearity  of  the  detector  was  determined  by  exposing  the  detector  to  a  wide  range  of  uniform 
x-ray  exposures  for  each  of  the  radiographic  techniques  described  above.  The  average  pixel 
values  were  computed  from  a  14.3  x  14.3  cm  region  located  near  the  chest  wall  section  of  the 
detector.  From  this,  the  relationships  between  mean  pixel  value  and  exposure  were  ascertained. 

2.4  Modulation  Transfer  Function 

The  presampled  MTF  was  measured  using  an  edge  method  similar  to  that  reported  elsewhere.4"9 
A  sharp  edge  test  device,  consisting  of  a  polished  0.1  mm  platinum-iridium  edge,  was  placed  in 
contact  with  the  detector  at  1  cm  from  the  chest  wall  edge  of  the  detector.  The  device  was 
oriented  with  a  5-10  degree  angle  with  respect  to  the  pixel  array.  An  image  of  the  edge  device 


was  then  acquired  using  an  exposure  of  16.2  |xC/kg  (62.6  mR).  The  presampled  modulation 
transfer  function  (MTF)  was  then  computed  from  the  edge  image  using  a  method  described  in  a 
previous  publication.10  The  MTF  was  computed  along  two  orthogonal  directions — the  Chest 
Wall-Nipple  (CN)  axis  and  the  Left-Right  (LR),  as  shown  in  Figure  1. 


Chest  Wall-Nipple 

fr - 

Detector 


_ ^  Left- Right 

Fig.  1.  Coordinate  system  for  measurements 

2.5  Normalized  Noise  Power  Spectrum 

To  characterize  the  system  noise,  flat-field  images  were  acquired  by  exposing  the  detector  to  a 
uniform  x-ray  beam.  The  exposure  was  simultaneously  measured  with  the  ionization  chamber 
reported  above.  The  Normalized  Noise  Power  Spectrum  was  then  computed  from  the  flat-field 
images  using  previously  published  methods.9, 10 

2.6  Detective  Quantum  Efficiency 

The  Detective  Quantum  Efficiency  (DQE)  was  computed  using  the  following  equation: 


DQE{u)  = 


MTF2  (u) 
q^deal  ’  E  '  NNPS{u) 


(1) 


where  MTF(u)  represented  the  presampled  modulation  transfer  function  measured  above,  qued 
corresponded  to  the  signal  to  noise  (SNR)  ratio  per  unit  exposure  for  an  ideal  energy-integrating 
detector,  E  was  the  exposure  value  at  the  detector  face  at  which  the  Normalized  Noise  Power 
Spectrum,  NNPS(u),  was  measured.11, 12  The  queai  was  computed  with  an  x-ray  simulation 
program  (xSpect,  Henry  Ford  Health  System)  that  utilized  a  semiempirical  model  to  simulate  the 
x-ray  spectra13  and  attenuation  effects.5 

2.7  Image  Lag  Measurement 

The  magnitude  of  multiplicative  lag  was  characterized  using  the  procedure  described  in  IEC 
62220-1. 14  First,  an  image  was  acquired  of  a  uniform  radiation  field.  The  second  image  was 
then  acquired  of  an  edge  device.  After  a  specified  delay  time  At,  a  third  image  was  acquired  of  a 
uniform  radiation  field.  The  image  data  were  then  examined  for  two  regions  within  the  images. 
The  first  ROI,  ROIi,  was  placed  in  an  area  of  the  images  that  did  not  contain  the  edge  device  in 
image  2.  The  second  ROI,  ROI2,  was  placed  in  an  area  that  was  inside  the  region  covered  by  the 
edge  device  in  image  2. 


Time  1 


Time  2 


Time  3 


ROI  1 


ROI  2 


Fig.  2.  Description  of  the  lag  measurement  procedure. 


Edge  Device 


The  detector  was  judged  to  have  acceptable  lag  effects  for  time  delay  At  if  it  passed  the 
following  criterion14 


(ImagelR0Ii  - ImagelROi2) - (Image3ROI1  - Image3ROi2)|  <  ^ 
Image  1R0I2  +Image3ROI2 
~2 


(2) 


3.  Results 


Fig.  3.  Plot  of  mean  pixel  value 
versus  exposure.  The  system 
showed  a  very  linear  response  with 
r2>  0.999. 


The  large  area  transfer  characteristics  of  the  detector  are  shown  in  figure  3.  The  detector 
maintains  its  linearity  over  two  orders  of  magnitude  in  exposure.  The  MTF  is  shown  in  figure  4 
for  the  CN  and  LR  axis.  The  MTF  along  these  axes  diverge  at  higher  spatial  frequencies.  The 


MTF  curves  are  summarized  in  Table  I  for  each  axis. 
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and  LR  axes 


RQA-M2 

(CN  Axis) 

RQA-M2 

(LR  Axis) 

0.2  MTF 

11.1  mm'1 

12.7  mm'1 

0.1  MTF 

12.8  mm'1 

14.8  mm"1 

0.5  mm'1 

0.983 

0.991 

2.5  mm'1 

0.858 

0.877 

5.0  mm'1 

0.65 

0.689 

Table  I.  Summary  of  the  detector’s  MTF  properties 


The  radial  traces  of  the  NNPS  multiplied  by  exposure  are  shown  in  figure  5  for  each 
radiographic  technique.  The  product  of  NNPS  and  exposure  should  remain  constant  for  strictly 
quantum  noise-limited  detectors.  However,  the  system  exhibited  exposure  dependencies.  For 
lower  exposures,  the  magnitude  of  this  metric  first  decreased  and  then  increased  with  increasing 


exposure. 
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Fig.  5.  Radial  trace  of  NNPS  for  various  exposure  levels 


Figure  6  shows  the  measured  DQE.  The  DQE  curves  showed  a  decline  at  low  frequency,  which 
was  expected  from  the  strong  low-frequency  component  of  the  NNPS.  As  well,  the  DQE 
increased  with  exposure  for  lower  exposure  values,  reached  a  peak  value,  and  then  decreased  for 


higher  exposures.  This  was  also  expected  from  the  behavior  of  the  NNPS  with  exposure.  The 
DQE  is  summarized  in  table  II. 


Fig.  6  Plots  of  the  DQE  at  various  exposures  along  CN  axis,  left,  and  LR  axis,  right. 


CN  Axis 

LR  Axis 

0.15  mm1 

46% 

47% 

2.5  mm'1 

49% 

50% 

5.0  mm'1 

31% 

34% 

Peak 

55% 

1 .25  mm'  1 

53% 

1.45  mm'1 

The  results  from  lag  measurements  are  summarized  in  table  III.  The  images  were  acquired  in  the 
order  indicated  in  table  3,  with  shorter  delay  time  tests  preceding  longer  delay  time  tests.  In 
general,  the  image  lag  for  the  detector  passed  the  test  established  by  the  IEC  (Eq  2).  However, 
an  interesting  phenomenon  occurred  for  the  75  pGy  exposure  with  5  minute  delay.  A  200  |lGy 
exposure  was  acquired  10  minutes  before  this  exposure.  It  appeared  that  this  high  exposure  still 
affected  the  detector  after  10  minutes,  as  a  75  (iGy  exposure  should  not  have  caused  a  larger  lag 
contribution  after  a  5  minute  decay  time  than  it  would  after  a  3  minute  decay  time. 


Table  III.  Summary  of  Multiplicative  Lag  Measurements 

Exposure  (|xGy) 

Decay  Time  (min) 

Metric 

Acceptable? 

75 

3 

0.002 

Yes 

150 

3 

0.0048 

Yes 

200 

3 

Ten  Minute  Wait 

0.0044 

Yes 

75 

5 

0.047 

No 

150 

5 

0.013 

No 

200 

5 

0.0022 

Yes 

4.  Discussion 

This  prototype  detector  has  excellent  resolution  properties,  as  shown  by  its  MTF.  There 
appeared  to  be  an  asymmetry  in  the  MTF,  as  it  diverged  for  the  CN  and  LR  axes.  As  the  edge 
device  was  placed  directly  on  the  detector  surface,  it  appeared  unlikely  that  the  focal  spot  caused 
such  asymmetries.  Future  work  is  needed  to  understand  the  cause  of  this  asymmetry.  The 
prototype  showed  structured  noise  contributions,  which  led  to  a  strong  low-frequency 
contribution  to  the  NNPS.  This  structured  noise  also  affected  the  DQE,  in  that  the  DQE  curves 
had  a  peak  and  then  decreased  for  lower  frequencies.  Finally,  image  lag  appeared  to  be  within 
the  parameters  established  by  IEC  62220-1, 14  but  high  exposures  led  to  unusual  behavior  in 
signal  retention,  even  after  a  long  decay.  This  prototype  showed  excellent  promise  and  it  is 
expected  that  future  work  will  correct  the  observed  structured  noise  and  lag  phenomena  with  a 
more  robust  calibration  technique. 
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